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1.0  EXECUTIVE  SUMMARY 


DoD  has  a  need  for  very  sensitive  midwave  and  long  wave  Infrared  photodetectors.  It  has 
long  been  recognized  that,  theoretically,  above  IK  quantum  photodetectors  offer  superior 
performance  over  thermal  photodetectors  (bolometers).  Accordingly,  with  the  discovery  of  high 
temperature  superconductivity  (HTS)  efforts  were  started  to  use  these  new  materials  for  making 
sensitive  midwave  and  long  wave  Infrared  quantum  photodetectors.  These  efforts  reflect  the 
technological  difficulty  and  high  cost  of  making  IR  photodetectors  in  semiconductors. 
Consequently,  the  expectation  is  that  with  superconductor  materials  more  sensitive  photodetectors 
with  longer  cut-off  wavelength  will  become  feasible. 

This  program  was  started  with  the  goal  of  demonstrating  Long  Wave  Infrared  quantum 
photodetectors  made  in  the  newly  discovered  HTS  materials.  Given  the  relative  immaturity  of  the 
HTS  field,  we  encountered  many  different  theoretical  and  experimental  problems  and  had  to 
formulate  many  solutions  to  these. 

We  have  formulated  a  quantum  photodetector  based  on  the  kinetic  inductance  of  a 
superconductor  and  named  it  the  Quantum  Superconducting  Kinetic  Inductance  Photodetector 
(QSKIP).  The  performance  of  this  QSKIP  is  analyzed  in  section  2.2.  The  analysis  reveals  that  the 
QSKIP  photoresponse  will  be  limited  by  the  Cooper  pair  binding  energy  2A,  believed  to  be  about 
32|am,  2A  of  YBCO.  The  QSKIP  photoresponse  and  sensitivity  are  computed  from  the  minimum  of 
the  Hamiltonian  energy  functional  and  linearized  Rothwarf-Taylor  equations.  Photoresponse  and 
sensitivity  expressions  are  computed  in  terms  of  the  quasiparticle  lifetime  and  indicate  Backgroimd 
Limited  Infrared  Performance  (BLIP)  at  very  low  photon  flux  levels.  At  low  temperatures  and  under 
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BLIP  conditions,  the  photoresponse  is  proportional  to  the  number  of  absorbed  photons.  Operating 
the  QSKIP  below  Hd  and  0.5  Jc  limits  the  noise  sources  to  fluctuations  in  the  condensate 
population.  A  Noise  Equivalent  Power  estimate  for  a  YBCO  QSKIP  is  about  7x10'*^  Watts  at  9K. 
These  predictions  are  based  on  eharacteristie  lifetimes  for  quasiparticle  generation,  quasipartiele 
recombination,  anharmonic  phonon  decay  and  phonon  trapping  in  a  BCS  type  superconductor. 

The  QSKIP  fabrication  is  described  in  section  2.3.  The  experimental  approach  and 
measurements  for  the  QSKIP  are  presented  in  seetion  3.0.  To  minimize  development  and  fabrication 
risk,  we  decided  to  use  a  YBCO  photodeteetor  with  a  niobium  based  SQUID  readout  circuit.  A 
generic  experimental  set  up  was  developed  for  measuring  photoresponse  of  superconducting 
photodetectors  from  1pm- 100pm  over  a  very  wide  temperature  range  (5-300K)  and  this  is  presented 
in  section  3.0. 

We  have  measured  the  characteristics  of  the  Nb  SQUID,  which  operated  properly  with  a 
eritieal  eurrent  of  0.2mA  at  5K.  The  YBCO  photodetector  material  was  supereonducting  at  a 
reduced  temperature.  The  as  deposited  YBCO  had  a  transition  temperature  of  about  90K  and  this 
was  reduced  to  about  60K  after  the  addition  of  the  Niobium  SQUID  readout  circuits.  We  have 
investigated  the  cause  of  reduction  in  the  transition  temperature  and  eoneluded  that  silicon  dioxide 
insulators  are  leaehing  oxygen  and  poisoning  the  YBCO.  Reoxygenation  of  the  YBCO  did  not 
reeover  the  90K  transition  temperature.  It  is  eoneluded  that  the  quasiparticle  lifetime  in  the 
proeessed  YBCO  was  too  degraded  to  show  photoresponse.  We  are  in  the  process  of  replacing  the 
silicon  dioxide  insulator  with  epitaxial  SrTi03.  This  work  is  eontinuing  on  Westinghouse  IR&D 
with  the  expectation  that  the  new  QSKIP  structure  will  produce  the  projected  performance  of  about 
10*^’ Watts  NEP. 
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2.0  QUANTUM  SUPERCONDUCTING  KINETIC  INDUCTANCE  PHOTODETECTORS 


In  this  chapter  we  provide  a  theoretical  framework  for  superconducting  quantum 
photodetectors  and  introduce  a  photodetector  structure  suitable  for  making  superconducting 
quantum  photodetectors. 

2.1  BACKGROUND 

The  pursuit  and  the  development  of  photodetectors  is  a  very  active  area  of  research.  This 
research  covers  a  very  broad  spectral  band,  from  y-ray  all  the  way  up  to  very  long  wave  infrared. 
The  operating  conditions  and  performance  of  photodetectors  operating  in  these  bands  varies 
according  to  the  photodetector  materials,  photodetection  method,  and  the  operating  spectrum.  Of 
particular  interest  are  photodetectors  operating  in  the  midwave  3-5pm  (MWIR)  and  longwave  8- 
12pm  (LWIR)  infrared  spectra.  Interest  in  these  spectral  bands  is  large  because:  (1)  the  atmosphere 
is  transparent  in  these  spectral  bands,  and  (2)  objects  at  300K  produce  significant  black  body  photon 
radiation  in  the  MWIR  and  LWIR  spectral  bands.  With  a  significant  black  body  photon  radiation, 
passive  night  imaging  and  target  discrimination  are  possible,  making  MWIR  and  LWIR  sensor  very 
important  for  DoD  and  commercial  applications. 

Detectors  fall  into  two  broad  categories:*’^  (1)  thermal  photodetectors  and  (2)  quantum 
photodetectors.  Thermal  photodetectors  absorb  the  incident  photon  flux  energy  with  an  optically 
“black”  coating  and  hence  are  sensitive  over  a  very  broad  spectral  band.  Thermal  photodetectors  are 
equilibrium  photodetectors  where  the  lattice  and  the  electrons  are  in  thermal  equilibrium.  The 
properties  of  the  photodetector  are  very  temperature  dependent  and  photoabsorption  produces  a 
temperature  change  that  generates  a  signal.  Quantum  photodetectors  require  an  energy  quanta  to 
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produce  an  excitaion  and  the  value  of  this  energy  quanta  translates  to  photoresponse  over  a  limited 
spectral  band.  Quantum  photodetectors  are  nonequilibrium  photodetectors  where  the  electrons  are 
not  in  thermal  equilibrium  with  the  lattice.  The  photodetector  properties  are  not  very  temperature 
dependent  and  photoabsorption  produces  quantum  excitations  in  the  electronic  states  and  thus  a 
photosignal.  Theoretically,  above  IK,  quantum  photodetectors  are  more  sensitive  than  thermal 
photodetectors  because  the  noise  level  in  quantum  photodetectors  is  smaller.  The  noise  in  quantum 
photodetectors  corresponds  to  Poisson  fluctuations  in  the  number  of  absorbed  photons,  whereas  the 
noise  in  thermal  photodetectors  corresponds  to  fluctuations  in  the  absorbed  energy.  This  basic 
difference  in  noise  mechanism  can  make  quantum  photodetectors  potentially  10^  more  sensitive  than 
thermal  photodetectors. 

The  pursuit  of  better  quantum  photodetectors  has  been  an  ongoing  endeavor  and  the 
discovery  of  High  Temperature  Superconductivity  (HTS)  widens  the  possibilities.  Just  before  the 
discovery  of  HTS,  Enomoto  and  Murakami^  made  photoresponse  measurements  on  granular 
BaPbo  7610  303,  and  reported  encouraging  results.  With  the  discovery  of  HTS  by  Bednorz  and 
Muller"*  in  1986,  many  proposals  were  submitted  for  applying  this  discovery  to  MWIR  and  LWIR 
photodetectors  operating  at  liquid  nitrogen  temperatures.  More  experimental  results  ’  ’  followed 

D 

reporting  on  the  photoresponse  of  various  HTS  photodetector  structures.  Forrester  et  al.  showed 
that  the  reported  photoresponse  signals  were  consistent  with  a  bolometric  response,  much  less 
sensitive  than  quantum  IR  photodetectors.  However,  because  of  technological  limitations  in 
existing  semiconductor  based  quantum  photodetectors,  in  particular  beyond  12|xm,  the  search  for  a 
superior  superconducting  photodetector  persists.^’ The  cited  approaches  for  realizing 
superconducting  photodetectors  have  limitations  that  include  poor  area  efficiency”  (or  quantum 


PARTI  II 


4 


efficiency),  complicated  readout  circuits,*^  and  speculative  photodetection  mechanism.^  To 
overcome  these  limitations,  a  simple  photodetector  structure  and  readout  circuit  are  needed.  These 
issues  were  addressed  in  this  program  and  are  reported  on  in  this  final  report. 

2.2  THEORY  OF  OPERATION 

We  describe  a  photodetector  structure  which  overcomes  the  aforementioned  limitations  and 
is  projected  to  acheive  Background  Limited  Infrared  Performance  (BLIP)  for  very  low  background 
radiation,  with  an  NEP  of  less  than  10’^^  Watts.  The  quantum  photodetector’s  photoresponse  is 
based  on  the  Cooper  pairs’  kinetic  inductance  and  hence  the  label:  Quantum  Superconducting 
Kinetic  Inductance  Photodetector*^  (QSKIP).  This  photodetector’s  photoresponse  is  sensitive  to 

photoinduced  changes  in  the  total  Cooper  pair  population  and  is  unlike  Josephson  junction  based 
photodetectors  that  are  only  sensitive  to  changes  in  the  condensate  population  within  a  coherence 
distance  ^  from  the  Josephson  Junction. 

The  QSKIP  is  intended  to  operate  at:  (1)  low  temperatures,  (2)  low  background  radiation,  (3) 
below  Hci,  (4)  below  0.5  Jc,  (5)  under  static  nonequilibrium  conditions,  and  (6)  not  under  pulsed 
laser  radiation,  as  reported  in  many  experiments,  ’  ’  but  modulated  at  low  frequencies  (<100KHz). 
The  reasons  for  choosing  such  operating  conditions  are  several.  In  quantum  photodetectors,  the 
photoresponse  increases  with  the  effective  quasiparticle  lifetime.  The  effective  quasiparticle  lifetime 
is  maximized  by  operating  at  low  temperature  and  near  thermal  equilibrium.  The  lowest 
photodetector  noise  is  achieved  by  operating  below  Hd  and  at  low  temperatures  (T/T^  <  0.1  ) 
because  noise  from  the  formation  and  movement  of  fluxoids  is  severely  inhibited.  Finally,  we 
predict  excellent  performance  for  the  QSKIP  when  operating  in  a  low  background  environment,  i.e. 
detection  of  cold  objects  (<  200K)  against  a  cold  backgrovmd  (<  200K) 
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The  QSKIP  is  made  from  a  simple  superconducting  film  patterned  into  a  closed  loop  with 
input  and  output  current  leads.  To  facilitate  the  photoresponse  and  sensitivity  calculation,  the 
QSKIP  is  represented  by  a  circular  loop*^’*^  illustrated  in  Fig.  2.1.  For  mavirmim  coupling 
efficiency  and  sensitivity,  the  readout  SQUID  is  directly  coupled*^  with  the  QSKIP,  see  Figure  2.2. 
In  operation,  the  QSKIP  is  biased  by  a  constant  dc  current  Iq,  Iq  <  0.5  Iq.  Current  Iq  is  divided  by  the 
QSKIP  circuit  into  currents  Ij  and  I2.  With  photoabsorption,  the  division  of  Iq  into  Ij  and  I2  changes 
and  this  change  represents  the  QSKIP  photoresponse. 


(b) 


Figure  2.1.  A  superconducting  loop  (a),  thinner  than  the  London  penetration  depth,  is  subjected  to  a 
bias  current  Iq  injected  at  point  “A”  and  removed  at  point  “B”.  The  bias  current  Iq  is  divided  into 
currents  Ij  and  I2  flowing  in  branches  1  and  2,  respectively.  The  QSKIP ’s  equivalent  circuit  (b) 
contains  the  magnetic  (Li=  Ln  -  Lj2,  L2=L22  -  L]2)  and  kinetic  (^1,^2)  inductance  terms  for  branches 
1  and  2,  respectively.  The  self  and  mutual  magnetic  inductances  are,  respectively,  represented  by 
(Lii,  L22)  and  (L12  =  L21).  li  is  represented  as  a  variable  inductor  to  indicate  changes  caused  by 
photoabsorption  in  branch  1 . 

The  current  division  in  the  QSKIP  is  calculated  from  the  minimum  of  the  superconductor’s 
energy  functional,  the  Hamiltonian.  The  condensate’s  Hamiltonian  is  expressed  in  terms  of  the 
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QSKIP  operating  currents,  Ij  and  I2,  and  the  Cooper  pair  density  Pcp(X)  is  estimated  to  be  constants 
Pcpi  and  Pcp2  in  branches  1  and  2,  respectively.  After  incorporating  the  vector  potential  A(X)  into 
the  Hamiltonian,  we  calculate  from  the  energy  extrema  the  QSKIP’ s  equilibrium  (or  static 
nonequilibrium)  operating  conditions  that  are  characterized  by  currents  and  I2.  Photoabsorption  of 
a  constant  photon  flux  in  only  one  branch  reduces  the  number  of  Cooper  pairs  in  that  branch 
resulting  in  a  new  static  nonequilibrium  state.  With  each  new  level  of  photoillumination,  the 
QSKIP ’s  current  division  changes  in  response  to  the  new  static  nonequilibrium  conditions.  The 
QSKIP’s  performance,  expressed  in  terms  of  the  Noise  Equivalent  Power^^  (NEP),  is  calculated  by 
combining  expressions  for  the  QSKIP’s  photoresponse  and  noise. 


OSKID  bias  current 


(a) 


■1A 


Figure  2.2.  The  YBCO  QSKIP^^  (a)  is  shaped  into  a  serpentine  pattern  to  increase  it’s  kinetic 
inductance  and  improve  isolation  againts  quasiparticle  diffusion  between  the  two  branches.  The 
readout  SQUID  is  directly  coupled  to  the  QSKIP  at  the  bottom.  The  photon  flux  only  illuminates 
the  right  half  of  the  QSKIP.  TTie  equivalent  circuit  (b)  for  the  symmetrical  QSKIP  and  readout 
SQUID  represents  the  kinetic  and  magnetic  inductances  for  branches  1  and  2.  The  QSKIP  is  biased 
with  a  dc  current  I^  and  the  SQUID  is  biased  with  a  dc  current  igQ. 
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2.2.1  THE  QSKIP  ENERGY  MINIMUM 


The  QSKIP  ring ,  illustrated  in  Fig.  2.1,  receives  a  dc  bias  current  Iq,  at  point  A,  and  divides 
it  into  currents  Ij  and  I2,  in  branches  one  and  two,  respectively,  and  Iq  exits  at  point  B.  The  total 
energy  of  the  condensate  Egg^ip  is  computed  by  integrating  the  energy  of  each  pair  over  the  pair 
density  Pcp(X)  and  the  QSKIP  volume  2V.  The  value  of  currents  Ij  and  I2  satisfy  the  energy  extrema 
of  the  superconductor’s  Hamiltonian  operating  on  the  condensate,  represented  by  a  wave  function 
'P.  The  Hamiltonian  written  in  terms  of  canonical  momentum*^  (the  squared  term  in  Eq.  1)  and  the 
Cooper  pairs  interacting  potential  (last  term)  is: 

»=^[p(X)+qepA(X)J  +  qe,<I)(X)+q„V„(X)  (1) 

Where  p(X)=nicpV  is  the  kinetic  momentum,  A(X)  is  the  vector  potential,  Vcp(X)  is  the  attractive 
potential  for  a  Cooper  pair  at  X,  0(X)  is  the  electrostatic  potential  energy,  qcp  is  the  electric  charge 

17 

for  a  pair  and  nicp  is  the  mass  of  a  pair.  We  simplify  equation  1  by  making  use  of  the  Coulomb 
gauge  V»A(X)=0,  and  the  so  called  London  gauge  AJ(X)=-A(X),  where  A=m<;;p/[pcp(X)qcp  ]  and 
J(X)  is  the  current  in  the  superconductor. 

The  Hamiltonian  in  Eq.  1  is  simplified  to  a  form  similar  to  the  BCS  Hamiltonian^^  i.e., 


1 

2nicp 


[p(X)]=-t 


— —  A(X)- J(X)  +q„®(X)+q„V„(X) 

^Pcp 


(2) 


We  have  used  the  expression  J(X)=qcppcp(X)p(X)/mcp  to  simplify  equation  1  and  obtain  an 
expression  for  the  Hamiltonian  with  two  terms  in  addition  to  the  BCS  terms:  the  second  term 
accounts  for  the  energy  present  in  a  flowing  superconducting  current,  while  the  last  term  with  0(X) 


PARTI  II 


8 


was  eliminated  by  BCS*^  by  solving  the  problem  relative  to  the  Fermi  Energy,  Ep.  Unlike  the  first 
two  terms,  the  last  two  terms  in  Eq.  2  are  not  a  function  of  currents  Ij  or  and  they  will  drop  out 
when  the  energy  minimum  of  the  QSKIP  is  calculated  relative  to  Ij  and  I2.  We  include  them  for 
generality. 

The  QSKIP  energy  Eqskip  is  obtained  by  integrating  over  the  Cooper  pair  density 
Pcp(X)='P*T'  where  this  density  spans  the  domain  of  the  Hermitian  operator  Incorporating 
pcp(X)='P*'P  and  substituting  J(X)=qcpPcp(X)p(X)/mcp»  fhe  Hamiltonian  given  by  Eq.  2,  the 
expression  for  the  QSKIP  energy  Eqskip  becomes: 

Eqskip=  J'P‘??^d^x  =  i  jAJ(X).J(X)d'x+i  jA(X)-J(X)d^x  + 

2V  ^  2W  ^  2V 

(3) 

jpcp(X)[<D(X)  +  Vcp(X)]d^x 

2V 

The  QSKIP  energy  Eqskip  Equation  3  has  three  terms:  (1)  the  first  term  represents  the 
condensate’s  kinetic  energy,  (2)  the  second  term  represents  the  energy  stored  in  the  field  produced 
by  a  flowing  supercurrent,  and  (3)  the  third  term  represents  the  potential  energy  terms. 

Since  the  Eqskip  energy  minimum  is  calculated  in  terms  of  currents  Ij  and  I2  [or  the 
corresponding  current  densities  Ji(X)  and  J2(X)]  we  need  to  express  Equation  3  in  terms  of  these. 
This  is  facilitated  by  doing  the  volume  integration  over  branches  1  and  2.  The  first  and  last  terms  in 
Equation  3  will  each  yield  an  additional  term  while  the  middle  term  will  yield  three  additional 
terms.  The  vector  potential  term  expressed  in  terms  of  volume  integral  over  branches  1  and  2  is: 
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Incorporating  these  substitutions  into  Eq.  3,  we  obtain  an  equation  for  the  QSKIP  energy  in 


terms  of  Ij  and  I2  and  the  kinetic  Zj  and  magnetic  Ly  inductance,  where  the  subscripts  “i”  and  “j” 
refer  to  branches  one  or  two,  i.e., 

EqSKID  ^“^2  ^1  2  ^ 


+  Jpc,(X)[0(X)+Vep(X)]d>x 


Where  =  (nicp/q  cp)Oi/(SiPcpi))  [4=(nk:p/q  cp)(y(S2PcP2))]  is  the  kinetic  inductance  of  branch 
number  one  (two)  with  a  Cooper  pair  density  pcpi  [Pcp2]»  a  cross  sectional  area  Sj  (S2)  and  a  length 
1,  (I2).  The  inductance  of  branch  one  (two)  is  labeled  by  Ln  (L22)  and  the  mutual  inductance 
between  the  branches  is  labeled  as  Li2=  L21,  by  symmetry. 

The  minimum  of  Egg^ip  is  calculated  by  taking  the  differential  of  equation  5  with  respect  to 
Ij  and  I2  to  yield, 


•=[A  +  “  [4  +  4]4  -  0 


We  make  use  of  the  fact  that  the  currents  flowing  in  the  QSKIP  branches  divide  Io=Ii+l2 
hence  5/5Ii=(5l2/5Ii)(l/5l2)=  -5/512  and  the  magnetic  inductance  terms  are  combined  into  effective 
inductance  terms;  Li=Lii-Li2  and  L2=L22-Li2.  It  should  be  noted  that  for  a  symmetrical  QSKIP 
Lii=L22  and  L2i=Li2.  Equation  6  holds  for  any  arbitrary  value  of  Cooper  pair  density  occurring  for 
each  level  of  constant  photodepairing. 
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2.2.2  QSKIP  PHOTORESPONSE 


In  this  section  we  calculate  the  QSKIP  photocurrent  from  the  constraints  governing  the 
photodetectors  operation  (given  by  Equation  6).  The  photocurrent  Ig  is  calculated  first  in  terms  of 
the  change  in  kinetic  inductance  Aij  (section  2.2.2.1).  Next,  using  the  modified  Rothwarf-Taylor 
Equations,  we  calculate  the  change  of  Mi  in  terms  of  changes  in  the  quasiparticle  density  ANq 
produced  by  photoabsorption  (section  2.2.2.2).  Finally,  we  relate  the  change  in  Mi  to  the 
quasiparticle  effective  lifetime  described  in  section  2.2.2.3. 

2.2.2.1  Photocurrent  Calculations 

The  photoresponse  signal,  Is=Ii-l2,  is  produced  as  the  division  of  Iq  changes  from  Ij  =  I2  to 
I,  I2  as  one  of  the  QSKIP  branches  is  exposed  to  a  constant  level  of  illumination,  and  the  currents 
change  to  satisfy  the  required  minimum  energy  conditions.  Without  loss  of  generality,  and  for  the 
sake  of  simplicity,  we  assume  a  totally  symmetric  QSKIP  geometry,  where,  under  no  illumination 
and  under  thermal  equilibrium,  the  Cooper  pair  density  Pcp(X)  is  the  same  in  branches  1  and  2 
producing  a  symmetric  current  division  i.e.  Ii=l2=Io/2.  Hence,  the  thermal  equilibrium  signal  current 
Is”,  and  the  associated  magnetic  flux  Ob”  inside  the  ring  are  zero.  Photoabsorption  in  branch  number 
1  only  breaks  the  QSKIP’s  symmetry  in  Pcp(X)  and  produces  a  photosignal  Is=  Ipla  and  a 
corresponding  net  magnetic  flux  Ob  inside  the  QSKJP’s  ring,  where  Ob=  Is(L,+L2).  Read-out  of  the 
photoinduced  signal  current  Ig  (for  directly  coupled)  or  Ob  (for  magnetically  coupled)  is  by  a 
SQUID. 

The  signal  current  Ig  is  produced  by  kinetic  inductance  changes  Mi  caused  by 
photoabsorption,  in  branch  1  only,  of  N<i,  photons  with  an  average  value  of  N^,”.  The  signal  current 
consists  of  a  circulating  current  Ig,  causing  a  corresponding  magnetic  flux  Ob.  The  signal 
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photocurrent  Is  directly  depends  on  Al^  and  is  computed  by  taking  the  differential  of  Eq.  6  to 


obtain: 


Is=—M,= 

A  ^  n  1 


"b  -^1  +0^2  -^2  ] 


-Ai, 


(7) 


Equation  7  is  valid  under  static  nonequilibrium  conditions  and  when  the  quasiparticle 
density  Nq  is  much  less  than  Pcp(X).  Using  partial  derivatives,  the  change  in  the  kinetic  inductance 
A£i  is  expressed  as: 


a/, 

apep(X)  dN 


Explicit  expressions  for  the  partial  derivative  factors  in  Eq.  8  are  obtained  using  Eq.  6,  the 
definition  of  (given  below  Eq.  5),  and  the  fact  that  each  Cooper  pair  forms  two  quasiparticles 
(5pcp(X)/oNq=-1/2).  To  complete  the  calculation  for  Is,  an  explicit  expression  for  ANq  is  needed 
that  represents  the  photoinduced  increase  in  the  quasiparticle  population  in  branch  1 . 


The  total  static  nonequilibrium  quasiparticle  density  Nq  is  expressed  as  a  sum  of  the  thermal 
equilibrium  quasiparticle  density  and  the  photoinduced  increase  in  the  quasiparticle  density, 
Nq=ANq+  Nq^^.  The  thermal  equilibrium  quasiparticle  density  in  each  branch,  based  on  a  BCS 
superconductor  at  T<0.5  Tc  ,  is  given  by:^* 


iVf  =  2iV(0)^|A(r)^,7’exp 


^  Ajjf 
V  k,T  j 


(9) 


where  N(0)  is  the  density  of  states  in  volume  and  energy  at  the  Fermi  surface,  ke  is  Boltzmann’s 
constant,  T  is  the  operating  temperature  in  degrees  Kelvin,  and  A(T)  is  the  temperature  dependent 
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superconducting  energy  gap.  An  expression  for  the  change  in  the  quasiparticle  density  ANq  from  the 
thermal  equilibrium  value  is  obtained  with  a  modified  form  of  the  Rothwarf-Taylor^^ 

equations. 


2.2.1.2  Photoinduced  Quasiparticle  Population 

The  Rothwarf-Taylor  equations,  10  and  11,  describe  a  superconductor  in  nonequilibrium  in 
terms  of:  the  quasiparticle  density,  Nq,  the  pair  breaking  (PB)  phonon  density,  and  the  absorbed 
photon  flux  density N,i„  as: 

dN„  2N„ 

— ^=2m,  riN,  +  — -  RN"  (10) 


dN,  ^  RN^  N, 

dt  2  Tg 


— +— 1[N„-N^] 

_^ES  "^AhJ 


2m,riN, 


(11) 


where  R=1/[tr^*^  Nq^^],  1/tr^^  is  the  thermal  equilibrium  quasiparticle  recombination  rate,  I/tr  is 
the  quasiparticle  generation  rate  by  PB  phonons  (since  only  PB  phonons  have  sufficient  energy 
2A(T)  to  break  Cooper  pairs).  The  phonon  depletion  occurs  by:  (1)  PB  phonon  escape  into  the 
substrate  with  rate  1/xes»  or  (2)  anharmonic  rate  decay  of  PB  phonons  into  multiple  phonons  (each 
with  energy  less  than  2A(T))  with  rate  1/tah. 


We  have  incorporated  into  the  Rothwarf-Taylor  equation  the  effects  of  photoabsorption  of  a 
photon  density  N<i„  acting  with  an  internal  quantum  efficiency  r|.  Specifically,  each  photoabsorbed 
photon,  with  energy  hv,  is  thermalized  by  electron-electron  and  electron-phonon  interactions. 
Quasiparticles  directly  photogenerated  by  electron-electron  interactions^'*’^^  are  represented  by  the 
term  2miqN(i„  in  Eq.  10.  Quasiparticles  generated  by  electron-phonon  interactions  are  represented 
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by  m2TiN0,  in  Eq.  11.  The  constants  ni]  and  m2  represent,  respectively,  the  number  of  quasiparticles 
and  PB  phonons  produced  by  photoabsorbed  photon  ’  .  In  these  calculations  we  assume  photons 
with  energy  (hv)  much  larger  than  2A  and  the  PB  phonons  with  energies  (hQ^A)  which  are  related 
by  an  approximation  hv«mi(2A)+m2(hQ).  Here  the  Fano  factor  has  been  neglected  .  Also  for 
YBCO  Tes  These  equations  are  simplified  into  linear  equations  by  invoking  static 

nonequilibrium  operating  condition  detailed  below. 

First,  since  the  QSKIP  is  operating  in  static  nonequilibrium  conditions  on  average 
dNQ/dt=dNQ/dt=dN(i/dt=0.  There  are  time  fluctuations  in  the  densities  of  N<i„  NQ,and  Nq  and  they 
manifest  themselves  as  noise  in  the  QSKIP,  and  these  are  considered  in  Section  2.2.3.  Here,  we  are 
only  interested  in  the  average  values  of  N,j>,  Nq,  and  Nq,  and  this  simplifies  Equations  10  and  1 1 
into  two  nonlinear  equations. 

Second,  the  nonlinearity  in  the  Rothwarf-Taylor  equations  is  also  removed  by  incorporating 
the  static  nonequilibrium  conditions  into  the  term  R  in  equations  10  and  11.  The  nonlinear  term 
represents  the  quasiparticle  recombination  rate  and  is  expressed  as  a  product  of  Nq  and  R.  The  Nq 
term  represents  the  number  of  ways  a  pair  can  be  formed  from  Nq  quaisparticles.  Specifically,  for 
Nq  quasiparticles  there  are  Nq(Nq-1)/2  ways  of  pairing  and  since  each  pairing  removes  two 
quasiparticles  we  obtain  the  Nq  factor,  for  Nq»1  .  The  term  R  represents  the  recombination  rate  for 
two  particular  quaisparticles.  The  term  R  is  not  equal  to  the  reciprocal  of  the  quasiparticle  lifetime 
1/tr^^,  calculated  by  Kaplan  et  al.^°  Kaplan  calculated  the  recombination  rate,  at  thermal 
equilibrium,  of  a  particular  quasiparticle  with  any  one  of  the  Nq^*^  available  quasiparticles.  Thus  at 
thermal  equilibrium,  the  factor  R=1/[tr^^q^^]  and  Nq^^  is  included  in  the  denominator  to  preclude 
double  coimting.  Rothwarf  and  Taylor  expressed  R  in  terms  of  the  thermal  equilibrium  values 
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and  Nq^^,  where  can  be  taken  from  Kaplan.^®  Since  the  QSKJP  operates  under  static 
nonequilibrium  conditions,  it  is  appropriate  to  express  R  in  terms  of  the  corresponding  static 
nonequilibrium  values,  specifically  R=1/[trNq].  This  representation  for  R  has  the  benefit  of 
linearizing  the  Rothwarf-Taylor  equations. 

The  solution  for  ANq  is  obtained  by  incorporating  the  above  conditions  into  equations  1 0 
and  1 1  and  solving  for  to  obtain: 


+  (12) 

TaH  "fTES 

This  indicates  that  photoabsorption  increases  the  PB  phonon  population  beyond  the  thermal 
equilibrium  number  and  the  increase  is  directly  proportional  to:  the  effective  internal  quantum 
efficiency  (mj+m2)q,  the  PB  removal  lifetime  (Trs  'tAH)/(T^ES'*^AH)»  the  photon  density  N^. 
Substituting  Eq.  12  into  Eq.  10  and  simplifying  by  recognizing  that  Tb«(xes  '^ah)/(''^es''^ah)»  we 
obtain 


N„'^  +  ANQ=2Ti(m,+  ITb) 


A 

Tah  Tes 
VTah  TTes/ 


Tr 


Xr 


N*  +2^N„ 

Tb  Xb 


EQ 


(13) 


In  Eq.  13,  Nq  is  replaced  by  Nq^^+ANq  to  provide  an  expression  for  the  change  in 
quasiparticles’  density  in  terms  of  and  Equation  13  is  valid  only  for  nonequilibrium 

conditions  where  Nq>  or  ANq>0.  This  requirement  is  satisfied  by  replacing  (in  Eq.  13) 
2(xR/xB)Nf2^^  by  Nq^^,  since  under  thermal  equilibrium  conditions,  when  Nq>=0,  the  nonequilibrium 
excess  quasiparticle  density  ANq=0.  Solving  Eq.  13  for  ANq  and  substituting  this  into  Eq.  8  and 
combining  with  Eq.  7,  the  explicit  expression  for  becomes 
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Is 


(14) 


Iq  I 

4p,,(X)(L,+i,) 

The  expression  in  the  square  brackets  in  Eq.  14  is  for  ANq.  The  photocurrent  Ig  depends  on 
the  effective  quantum  efficiency  [(mi+m2)T|],  the  effective  quasiparticle  lifetime 
'^r'^ah'*^es/[('^ah‘’^es)'^b]5  N<i,.  The  effective  internal  quantum  efficiency  term  is  made  up  of;  (1) 

the  internal  photon  absorption  efficiency  t]  within  the  superconductor,  (2)  the  number  of 
quasiparticles  generated  by  electron-electron  interaction  (mi),and  (3)  the  number  of  quasiparticles 
generated  by  electron-phonon  interactions  (m2).  Thus  each  absorbed  photon,  with  energy  several 
times  larger  than  2A  generates  up  to  hv/A  quasiparticles.  As  expected,  the  effective  quasiparticle 
lifetime  has  a  direct  influence  on  determining  the  signal  Ig  and  it  is  considered  next. 

2.2.23  Effective  Quasiparticle  Lifetime 

We  estimate  the  quasiparticle  lifetime  assuming  a  BCS  type  superconductor,  based  on  “s” 
wave  pairing.  This  is  appropriate  for  low  temperature  materials  but  is  presently  unresolved  for  high 
temperature  materials,  where  the  nature  of  the  superconducting  gap  is  still  in  question.  It  should  be 
noted  that  calculations  taking  account  of  “d”  and  “s”  wave  pairing  have  been  made  by  Sergeev  and 
Reizer  with  results  similar  to  our  calculations.  With  this  caveat,  we  proceed  to  estimate  the 
quasiparticle  lifetime  using  a  BCS  model. 

The  intrinsic  quasiparticle  lifetime  at  or  near  thermal  equilibrium  is  very  sensitive  to 
temperature  and  for  a  BCS  superconductor  at  T<  O.ST^  the  quasiparticle  lifetime  is  given  by 
Kaplan  as^*’: 


^ ,  'J  Tah  Tes  Tr 

r\  mi  +  nhl - N* 

VTaH  "flES/  Tb 
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V(T)=- 


V"  2A, 


k,T 


0.1  To  exp 


k,T 


where,  for  low  temperature  superconductors,  a  representative  values  for  Xq  calculated  by  Kaplan^^ 
et.  al  is  approximately  10'**  seconds.  At  small  reduced  temperatures  (T<0.1Tc)  and  for  a  BCS 
superconductor  (Aq  »  1.72  ),  the  intrinsic  recombination  lifetime  (computed  with  Eq.  15) 

becomes  very  large  (trw  SOpsec)  and  is  much  larger  than  Xq.  Since  high  temperature 
superconductors  are  projected  to  have  an  even  larger  energy  gap,  even  longer  quasiparticle  lifetimes 
are  expected.  For  example,  in  YBCO  the  energy  gap,  obtained  from  specific  heat  measurements 
(Ao»2.5kBTc),  is  larger  than  the  BCS  gap  leading  to  a  very  long  projected  intrinsic  quasiparticle 
lifetime,  70msec. 

The  effective  quasiparticle  lifetime  depends  on  several  components;  the  intrinsic 
quasiparticle  lifetime,  x^;  the  phonon  escape  lifetime  into  the  substrate,  Xps;  the  anharmonic  phonon 
decay  lifetime,  x^h;  and  the  pair  breaking  lifetime,  Xg.  The  increase  in  the  effective  quasiparticle 
lifetime  is  evident  from  the  solution  to  the  Rothwarf-Taylor  equations  12  and  13  where  the 
combination  of  lifetimes  Xr,  Xg,  Xgs,  and  x^jj  results  in  an  effective  quasiparticle  lifetime  Xgpp  given 


Teff 
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The  effective  quasiparticle  lifetime  is  limited  by  the  depletion  rate  of  nonequilibrium  PB 
phonons  through:  (1)  escape  into  the  substrate  with  effective  rate  [(tr/tb  )tes]  ^  and  (2)  anharmonic 
phonon  decay  with  effective  rate  [(tr/xb  )xah]  At  lower  temperatures  and  near  thermal 
equilibrium,  these  rates  are  much  slower  because  (tr/xb  )»1  • 

The  ratio  of  Xr/xr  is  estimated  in  terms  of  the  quasiparticle  density  Nq  and  the  total  Cooper 
pair  density.  Using  the  golden  rule  for  time  dependent  perturbation  theory, the  quasiparticle 
recombination  rate  is  given  as  l/xR=27i/h|<f|»?^iNli>P  FNq,  where  is  the  interaction 
Hamiltonian,  F  contains  coherence  and  occupation  factors,*^  Nq  is  the  density  of  available 
recombination  sites  for  going  from  the  initial  state  |  i>  to  the  final  state  |  f>.  For  recombination 
conditions  where  Pcp(X)»Nq.  Similarly,  using  the  golden  rule,  the  pair  breaking  rate 
l/xB=27r/h  I  <f  I  »?^iN  I  i>  I  ^  rPcp(Ef)  where  the  value  of  |  <f  I  »?^in  I  i>  P  is  the  same  squared  matrix 
element  used  for  calculating  1/xr.  Quasiparticle  generation  rate,  under  conditions  where 
Pcp(X)»Nq,  is  proportional  to  Pcp(X)»  since  this  is  approximately  equal  to  the  number  of 
unoccupied  quasiparticle  states.  From  these  expressions,  we  estimate  for  the  ratio  of  the 
quasiparticle  generation  to  the  quasiparticle  recombination  lifetime  as  Xr/xb»Pcp(X)/Nq. 

The  effective  anharmonic  decay  and  phonon  escape  lifetimes  are  greater  than  the  intrinsic 
lifetime  values,  x^h  and  because,  at  low  temperatures,  the  nonequilibrium  photodeposited 

energy  resides  most  of  the  time  in  the  form  of  excited  quasiparticles  (since  Xr/xb»Pcp(X)/Nq.»1), 
and  the  photodeposited  energy  can  only  escape  when  it  exists  in  the  form  of  nonequilibrium 
phonons  and  not  quasiparticles.  Since  the  photodeposited  energy  exists  for  the  largest  fraction  of 
time,  Xr/(xr+Xb),  in  the  form  of  excited  quasiparticles,  there  is  proportionally  very  little  time  during 
which  the  nonequilibrium  PB  phonons  can  escape  into  the  substrate  or  undergo  anharmonic  decay. 
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Therefore  at  low  temperatures,  the  effective  lifetimes  x^h  aiid  Xes  become  very  large.  With  a  longer 
Tepp,  the  QSKIP  expected  performance  is  significantly  improved  and  this  is  evident  from  the  Noise 
Equivalent  Power  (NEP)  calculations  in  section  2.2.5.  The  noise  current  51  is  calculated  next. 


2.2.3  NOISE  LEVEL  IN  QSKIP 

The  QSKJP’s  noise  current  51  is  calculated  from  the  Root  Mean  Square  (RMS)  fluctuations 
in  the  signal  current  Is,  specifically,  (5I)^=  <Is^>  -  <Is>^  where  the  brackets  <  >  represent  the 
expected  value.  A  superconductor  operating  in  the  current  state,  in  principle,  has  two  major  noise 
sources  possible  :  (1)  fluctuations  in  the  Cooper  pair  density  and  (2)  fluxoid  noise.  The  QSKIP  is 
intended  to  operate  at  low  temperatures  and  at  lo  ^  0.5  Ic*  Under  these  conditions,  the  operation  of 
the  QSKIP  can  be  readily  limited  below  Hci  thereby  minimizing  the  effects  of  fluxoid  noise. 
Accordingly  we  calculate  the  noise  in  the  QSKIP  as  primarily  due  to  fluctuations  in  the  Cooper  pair 
density.  Expressions  for  the  noise  current  51  are  derived  in  terms  of  the  statistical  fluctuations  in  the 
density  of  Cooper  pairs  pcp(X)  (or  quasiparticle  density  Nq)  in  branches  1  and  2  (see  Fig.  2.1).  The 
QSKIP  noise  current,  51,  is  expressed  in  terms  of  the  quasiparticle  density  variances,  CTqi  ,  and,  CTq2  , 
in  branches  1  and  2,  respectively.  The  expression  for  the  total  noise  current  51,  is  calculated  as  is  the 
signal  current  Is  by  using  Eq.  14  and  adding,  in  quadrature,  contributions  from  branches  1  and  2. 
Making  use  of  Eq.  7  and  8,  the  equation  for  the  QSKIP’s  RMS  noise  current  is: 
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Cooper  pair  density  (pcp(X)»Nq)  hence  the  probability,  P,  for  breaking  a  Cooper  pair  is  very  small 
and  we  estimate  the  average  quasiparticle  density  by  Ppcp(X)«Nq^‘^.  Under  these  conditions,  it  can 
be  shown  that  the  quasiparticle  statistics  can  be  represented  by  a  Poisson  Probability  Density 
Function  (PDF).  For  a  Poisson  PDF,  the  variance  c!q2  is  equal  to  the  mean  quasiparticle  density  in 
branch  2,  aQ2^=(l  A^)Nq^^.  The  mean  quasiparticle  density  is  given  by  Equation  8. 

Because  of  photoabsorption  ctqiV  The  quasiparticle  density  in  branch  1,  Nq,  is 
expressed  in  terms  of  the  thermal  equilibrium  quasiparticle  density  and  the  additional  density 
of  static  nonequilibrium  quasiparticles,  ANq,  generated  by  the  absorption  of  N^.  This  corresponds  to 
a  nonequilibrium  state  that  is  different  from  thermal  models  where  the  increases  in  the  quasiparticle 
population  are  due  to  heating  of  the  superconductor.  Here,  the  quasiparticles’  nonequilibrium  state 
is  maintained  by  an  excess  static  nonequilibrium  of  PB  phonons.  The  total  PB  phonon  density  is 
made  up  of  contributions  due  to  thermal  equilibrium  PB  phonons  represented  by  Bose- 

Einstein  statistics,  and  contributions  due  to  photogenerated  nonequilibrium  PB  phonons.  The 
photogenerated  nonequilibrium  phonons  span  a  narrow  energy  range  because  photoexcited  hot 
quasiparticles  thermalize  most  efficiently  by  the  emission  of  energetic  phonons,  to  which  they  are 
most  strongly  coupled.  These  energetic  phonons  are  primarily  the  PB  phonons  that  maintain  the 
nonequilibrium  quasiparticle  population  by  virtue  of  the  fact  that  they  participate  in  the  quasiparticle 
generation  and  recombination  process,  since  they  have  a  minimum  energy  of  «2A.  This 
nonequilibrium  phonon  representation  is  similar  to  (but  not  the  same  as)  the  T*  approach  used  by 
Parker.  Parker  represented  the  phonon  densities  in  a  superconductor  under  static  thermal 
nonequilibrium  as  a  sum  of:  (1)  thermal  equilibrium  phonon  density  represented  by  a  thermal 
equilibrium  temperature  Tq  and  (2)  an  excited  phonon  density  represented  by  a  higher  thermal 
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equilibrium  temperature  T*>To.  Here,  we  are  representing  the  static  nonequilibrium  phonon 
distribution  as  a  sum  of  two  phonon  densities:  (1)  the  thermal  equilibrium  PB  phonon  density, 
and  (2)  a  nonequilibrium  phonon  population  ANQ=Nn-  (see  Eq.  12),  which  span  a  narrow 
energy  range  with  an  energy  greater  than  2A.  It  should  be  clear  that  this  static  nonequilibrium 
distribution  of  quasiparticles  and  PB  phonons  cannot  be  described  by  two  thermal  equilibrium 
distributions  at  different  temperatures,  the  T  model. 


The  quasiparticle  variance  is  determined  by  the  processes  ongoing  in  branch  1.  In  static 
nonequilibrium,  the  effective  quasiparticle  lifetime  is  much  larger  than  the  intrinsic  quasiparticle 
lifetime  because  of  the  many  quasiparticle  generation  and  recombination  events  occurring  before  the 
photoabsorbed  energy  is  removed  with  the  anharmonic  decay  or  escape  of  the  PB  phonons.  These 
numerous  generation  and  recombination  events  (experienced  by  each  quasiparticle,  and  governed  by 
the  Roth warf-Tay lor  Eqs.  10  and  11)  effectively  decouple  the  quasiparticle  statistics  from  the 
photogenerating  photon  statistics.  Again  given  our  QSKIP  operating  constraints,  as  in  the  case  of 
the  unilluminated  branch  2,  the  quasiparticle  density  is  much  smaller  than  the  Cooper  pair  density 
(Pcp(X)»Nq).  This  corresponds  to  conditions  where  the  probability  for  breaking  a  Cooper  pair,  P, 
is  very  small,  thus  making  the  average  quasiparticle  densities  Nq^^+ANq«PNcp.  Under  these 
conditions,  the  quasiparticles  obey  Poisson  statistics  and  the  variance  agi  is  equal  to  the  mean 
quasiparticle  density  in  branch  1,  given  by  aQi^=(W)[NQ^^+ANQ]  (see  Eq.  13).  Using  the  explicit 
expression  for  the  quasiparticle  density  variances,  the  expression  for  the  QSKIP  noise  in  Eq.  17 
becomes: 


Ip  I,  |2ri(m,+m,)i:,ppN^+2N^*^ 

4Pcp(X)(/,  +L,)V  V 


(18) 
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The  square  root  term  includes  the  noise  from  the  ANq  additional  quasiparticles  photogenerated  in 
the  illuminated  branch  and  beyond  the  2Nq'^^  thermal  equilibrium  density  in  both  branches  . 


2.2.4  RESPONSIVITY  AND  NEP  OF  QSKIP 

The  QSKIP  operates  in  the  zero  resistance  current  state  and  accordingly  the  responsivity 
(91i)  and  NEP  (calculated  at  the  end  of  this  section)  are  calculated  in  terms  of  current  and  not 
voltage.  The  calculated  responsivity  is  readily  expressed  in  terms  of  the  photoinduced  magnetic  flux 
Ob  simply  by  multiplication  of  iRi  by  the  QSKIP  total  magnetic  inductance  L1+L2.  The  expressions 
for  the  current  responsivity  9Ii  is: 


31, 

Wattv 


1  Ic 


hv  vn; 


(19) 


The  explicit  expression  for  91i( A/Watt)  is  obtained  by  substituting  Is  (Eq.  14)  into  Eq.  19  and 


dividing  by  the  incident  photon  power,  VhoN,i,°,  to  yield: 


Vhv4p„(X)(e,  +L,) 


2Ti(m,  +  mj  )r 


EFF 


(20) 


The  responsivity  is  directly  proportional  to  the  effective  quasiparticle  lifetime  Teffj  as  is 
expected  for  a  photoconductive  nonequilibrium  quantum  photodetector.  At  low  temperatures,  the 
responsivity  is  greatly  increased  because  the  ratio  Tr/tb»Pcp(X)/Nq»1  (see  section  2.2.2.3  and 
calculations  by  S.  B.  Kaplan^”  et.  al)  resulting  in  an  effective  quasiparticle  lifetime  significantly 
larger  than  the  phonon  trapping  lifetime  Tes  (produced  by  acoustical  mismatch  at  the  film  substrate 
interface)  or  the  anharmonic  phonon  decay  lifetime  Tah-  Thus,  the  maximum  responsivity  occurs  at 
low  temperatures. 
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BLIP  is  achieved  when  Tgpp  is  sufficiently  large  so  that  the  photon  radiation  shot  noise  is 
equal  to  the  photodetector  noise  and  the  readout  circuit  noise,  added  in  quadrature.  From  Equation 
18,  it  is  evident  that  one  requirement  for  BLIP  is  that  Ti(mi+m2)TEppN^®>NQ^^.  The  NEP  is 
generically  represented  by: 

NEP(Watts)  = 

Under  BLIP  conditions,  the  photon  radiation  shot  noise  is  equal  to  the  photodetector’s 
internal  noise.  In  section  2.2.4  it  was  shown  that  the  photodetector’s  noise  can  be  represented  in 
terms  of  the  quasiparticle  variances  aqi  and  <Tq2  ,  that  obey  Poisson  statistics.  Hence,  noise  from 
the  photogenerated  quasiparticles  is  equal  to  the  noise  from  thermal  equilibrium  quasiparticles  when 
q(mi+m2)TEppN<i>®«NQ^^,  see  Eq.  18.  The  expression  for  the  QSKIP’s  BLIP  NEP,  obtained  with  Eq. 
20,  (and  assuming  negligible  electronic  readout  noise)  is 


NEPble, 


=  hv 


VNJ'^AF 

"b  ^2  EFF 


(22) 


The  QSKIP  NEP  varies  with:  VNq^^,  the  effective  internal  quantum  efficiency  [Ti(mi+m2)], 
the  electrical  operating  bandwidth  of  the  readout  circuit  AF,  and  the  effective  quasiparticle  lifetime, 
Tgpp.  The  VNq^^  term  corresponds  to  variance  of  the  thermal  equilibrium  quasiparticle  population 
(Poisson  statistics). 


2.2.5  The  NEP  of  a  YBCO  QSKIP 

The  NEP  is  computed  for  a  100pm  square  by  lOOnm  thick  YBCO  QSKIP.  In  computing  the 
NEP,  we  determine  the  operating  conditions  necessary  for  achieving  BLIP  performance. 
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Specifically,  three  conditions  are  required  and  are  satisfied:(l)  the  QSKIP  is  assumed  to  operate  at 
sufficiently  low  temperatures  such  that  the  photoinduced  quasiparticle  density  is  equal  to  or  larger 
than  the  thermal  equilibrium  quasiparticle  density  Nq^^;  (2)  the  QSKIP  bias  current  is  limited  to 
less  than  0.5Ic  to  insure  operation  below  Hd,  thereby  circumventing  noise  problems  associated  with 
fluxoid  generation  and  transport,  and  (3)  the  photon  density  levels  are  selected  so  that  the  noise 
associated  with  the  QSKIP  dominates  over  the  readout  noise,  i.e.  5I>5Isquid-  The  last  condition  is 
the  most  stringent  and  we  consider  it  first. 

To  satisfy  the  requirement  that  5I>5Isquid  several  operating  conditions  need  to  be  met:  (1) 
large  Teff>  (2)  lo  sufficiently  small  to  inhibit  fluxoid  noise,  and  (3)  Pcp(X)»Nq.  Presently  for 
YBCO,  there  are  no  universally  accepted  theoretical  or  experimental  values  available  for  Teff  and 
Pcp(X)  and  we  proceed  to  estimate  these. 

(1)  Estimating  Teff  in  a  YBCO  QSKIP  requires  knowledge  of  Tr/xb,  "Ces*  '^ah*  and  2A(0).  The 
value  of  2A(0)  is  needed  to  determine  the  PB  phonon’s  frequency  for  estimating  Tes>  and  Xah>  since 
an  estimate  already  exists  for  Xr/xb»Pcp(X)/Nq  (see  section  2.2.2,3).  From  specific  heat 
measurements^^  in  YBCO,  the  measured  pair  binding  energy  2A(0)  is  about  32meV»5kBTc.  This 
energy  gap  value  is  consistent  with  the  Raman  measurements  which  from  line  broadening  identified 
the  B,g  (out-of-phase  oxygen  plane  vibration)  at  340cm'*  as  being  particularly  strongly  coupled  with 
the  electrons'*'*’'**  and  this  is  inferred  as  being  related  to  superconductivity.  This  strongly  suggests 
that  the  PB  phonons  are  energetic  Bjg  optical  phonons  whose  energies  are  approximately  equal  to 
the  YBCO  Debye  energy'*^  phonons,  which  corresponds  to  the  value  of  the  energy  gap, 
2A(0)«32meV. 
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The  estimate  for  the  phonon  escape  lifetime  Xps  is  based  on  several  experimental  results. 
Experimental  measurements  with  acoustic  phonons  indicate  a  temperature  independent  99% 
scattering  probability  at  the  YBCO/substrate  boundary  (tes«10’  sec  for  lOOnm  YBCO  film)  .  For 
PB  phonons  there  are  no  similar  measurements,  however,  the  YBCO  phonon  dispersion  curves  in 
the  “c”  direction"^  indicate  that  there  are  no  phonon  propagating  modes  with  sufficient  energy  to 
break  Cooper  pairs.  Thus  from  a  c-axis  oriented  HTS  film,  PB  phonons  will  not  be  depleted  by 
escaping  into  the  substrate  because  of  the  relatively  poor  coupling  between  the  supercurrent 
carrying  CuO  planes.  Thus,  we  expect  that  the  depletion  of  PB  phonons  is  primarily  by  anharmonic 
phonon  decay,  i.e.  Tah«  '^es- 

Experimental  values  for  t^h  are  not  available  for  PB  phonons.  The  value  of  t^h  is  estimated 
from  measurements  on  YBCO  subjected  to  large  deviations  from  thermal  equilibrium,  and  induced 
by  excitations  with  energetic  laser  pulses^’.  At  large  deviations  from  thermal  equilibrium, 
Pcp(X)wNq  resulting  in  Tr/tb*!  (see  section  2.2.2.3),  and  the  measured  effective  relaxation  lifetimes 
should  correspond  to  TAH*10"100psec.  At  low  temperatures  and  near  thermal  equilibrium  thermal 
Xr/tb»Pcp(X)/Nq»10^,  leading  to  an  estimate  for  XEPF^lmsec. 

(2)  The  value  of  Iq  <  0.51^  is  chosen  to  facilitate  operation  of  the  QSKIP  below  Hd  and  outside 
the  fluxoid  formation  regime  and  associated  noise  effects.  A  reasonable  value  for  the  critical  current 
density  in  YBCO  at  low  temperatures  is  w  2  x  10^  A/cm^.  For  a  lOOnm  thick  YBCO  film,  7p,m 
wide,  IcwO.  14A  and  an  operating  current  Iq  <  0.1  Amp  is  less  than  Ic-  The  self  induced  magnetic  field 
produced  by  Iq  is  less  than  90  Gauss,  which  is  below  the  value  of  Hd  in  YBCO,'*^’'^^’'’’  thereby 
facilitating  the  maintenance  of  the  QSKIP  in  the  Meissner  state,  and  eliminating  flux  flow  noise. 
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(3)  The  Cooper  pair  density  for  YBCO,  Pcp(X),  is  estimated  from  specific  heat  measurement 
that  identify  the  contribution  of  the  electronic  component,  yp.  The  density  of  states  at  the  Fermi 
surface‘s*  N(0)w3Yp/(7r^kB^)  is  estimated  from  measured  values'*^  of  Yp  [between  1.4-2.9mJ/(g-atom- 
K^)]  to  be  about  2.8-5.8x10^’  states/(eV-m^).  Furthermore,  the  specific  heat  jump"*^  at  Tc  [AC/ypTc] 
varies  between  0.4-1.11  instead  of  the  BCS  value  of  1.43.  Loram  et.  al^^  measured  Yp  at  lOK  in 
YBCO  and  showed  it  vanishes,  hence  the  reduced  jump  in  the  specific  heat  at  Tc  is  not  attributed  to 
multiphase  YBCO  sample,  with  normal  regions,  but  it  is  simply  an  indication  of  a  smaller  fraction 
(10%-30%)  of  k-space  pairing  than  in  a  BCS  superconductor.  Accordingly,  our  estimate  of  Pcp(X) 
is  10%-30%  of  A(0)N(0)/2  «  0.2-1.4x10^^  pairsW.  For  a  serpentine  QSKIP  structure  with  7pm 
wide  lines  formed  in  a  0.1pm  thick  YBCO  film  and  within  a  50x100pm  area,  the  calculated  value 
for  V«5xl0'^^  m^,  which  corresponds  to  a  kinetic  inductance  of  about  .^i«1.4-8nH. 

We  can  now  proceed  to  calculate  the  value  of  N<i,®  needed  to  ensure  that  the  SQUID  read  out 
noise  is  less  than  the  QSKIP  noise,  5Isquid<  SI.  Good  modulation  of  the  readout  SQUID  output 
voltage  with  magnetic  flux  requires  LsqIsq  <  where  Igg  and  Lsq  are  the  SQUID  critical  current 
and  loop  inductance,  respectively.  Designing  a  SQUID  with  IsQ»100pA,  is  reasonable  given  a  loop 
inductance  LsQ»20pH.  These  SQUID  operating  conditions  yield  a  Y=hIsQ/ekBT»600  and  this  is 
consistent  with  operation  that  prevents  significant  noise  rounding  in  the  I-V  curves^®.  For  a  directly 
coupled  SQUID  readout  circuit,*^  and  under  BLIP  conditions,  we  require  that  SIsquid^  SI.  An 
achievable  SQUID  flux  noise  level^*’^^  is  «3xlO’’c>o/Hz’^«6xlO'^^  Webers/VUz,  and  for  LsQ«20pH, 
we  compute  5IsQuiD«3xl0’^*Amps/VHz.  To  satisfy  the  condition  SIsquid*^  SI,  we  need  to  select  a 
sufficiently  large  N,j)”  (see  Eq.  18)  to  satisfy  the  following  expression: 
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2ri(ni, +1112)1 


4p^^(X)(L.+i,)^ 


—  >3x10‘"VaF 


(23) 


Evaluating  Eq.  23  by  substituting  average  values  for:  2j«4.7nH,  ii/(Li+/,)»l,  Io»0.1  Amps, 
Pcp(X)w0.8xl0^^pairsW,  ri(mi+m2)«l,  and  30Hz  for  AF,  we  obtain  TEFFN<i,*^>1.4xlO**photonsW. 
For  a  projected  XEFF-lnisec,  we  determine  that  N<i,°>1.4xl0^'  photonsW-sec  (corresponding  to 
more  than  1.4x10*®  photons/sec-cm^)  is  required  for  photon  noise  to  dominate  over  the  readout 
noise,  i.e.,  61  >  6Isquid- 

For  BLIP  performance,  in  addition  to  the  requirement  that  51  >  5Isquid>  N,j,®  needs  to  be 
sufficiently  large  so  that  the  photodetector  is  dominated  by  the  photoinduced  quasiparticle  noise, 
that  is,  q(m,+m2)TEFFN(i,®>NQ®‘^  (see  Eq.  18).  Below  9K,  the  value  of  is  readily  calculated 
with  Eq.  9,  where  the  product  of  all  the  preexponential  factors  is  approximately  equal  to  about 
3.8x10^^  quasiparticlesW.  Evaluating  Eq.  9  at  9K,  with  A(0)«16meV,  we  obtain  Nq^^w1.6x10*® 
quasiparticlesW  and  this  leads  to  the  requirement  that  TeffN(i,”>  4x10*®  quasiparticles W 
[assuming  as  before  that  r|(mi+m2)»l].  This  is  a  less  strict  requirement  than  61  >  6Isquid  which 
requires  that  TeffN<d®>1  .4xl0**photons/m^  (see  bottom  of  previous  paragraph). 

For  6I>5Isquid»  ’^e  calculate  the  NEP  at  BLIP.  At  9K,  for  TEFF*lnisec, 
N<D®>1.4xlO^*photons/m^-sec  (corresponding  to  1.4x10*®  photons/sec-cm^).  For  O.leV  photons,  the 
equivalent  NEP  in  a  30  Hz  bandwidth  is  about  7x10’*’  Watts.  The  NEP  is  limited  by  the  readout 
circuit  noise.  Better  NEP  might  be  possible  with  smaller  operating  bandwidth,  or  larger  values  of 
Teff  or  A(T),  which  are  potentially  possible  with  HTS  material  with  a  higher  Tq.  The  QSKIP 
performance  is  currently  being  experimentally  investigated. 
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2.2.6  CONCLUSIONS  ON  QSKID  OPERATION 


We  have  described  a  Quantum  Superconducting  Kinetic  Inductance  Detector  suitable  for 
operation  in  a  low  background  environment  and  compatible  with  LTS  and  HTS  material 
technologies.  The  QSKIP  responsivity  and  NEP  are  analyzed  using  the  Hamiltonian  energy 
functional  and  linearized  Roth warf-Tay lor  equations.  The  QSKIP  is  intended  to  operate  at  low 
temperatures,  in  the  current  state,  and  below  Hd  to  maximize  photoresponsivhy  and  to  suppress 
flux  flow  and  Johnson  noise,  present  in  resistive  current  operation.  The  noise  remaining  is  due  to 
phonon  induced  fluctuations  in  the  Cooper  pair  population,  and  this  noise  is  significantly  reduced 
because  at  T«Tc  the  population  of  PB  phonons  (hQ  >  2A(T))  is  very  small.  As  in  semiconductor 
based  photoconductive  quantum  photodetectors,  the  QSKIP  responsivity  increases  with  the  effective 
quasiparticle  photoexcitation  lifetime.  Because  the  energy  gap  in  superconductors  (unlike  the 
bandgap  in  conventional  semiconductor  photodetectors)  is  comparable  to  the  energy  of  PB  phonons, 
trapping  of  these  phonons  increases  the  effective  quasiparticle  lifetime  and  thereby  maximizes  the 
QSKIP  photoresponsivity.  The  trapping  of  PB  phonons  in  YBCO  is  very  effective  because, 
according  to  experiment,  the  PB  phonons  cannot  travel  in  the  “c”  direction  (toward  the  substrate) 
since  such  large  energy  phonon  modes  do  not  exist  for  the  “c”  direction.  Thus,  the  PB  phonons  will 
be  effectively  prevented  from  escaping  into  the  substrate  and  can  only  escape  by  lateral  motion  or 
anharmonic  decay.  The  QSKIP  is  projected  to  have  a  spectral  photoresponse  up  to  2A(T),  the 
superconductor’s  energy  gap  «32p,m  for  YBCO  at  T«Tc.  Such  a  wide  spectral  response  is 
currently  achieved  only  with  thermal  photodetectors.  Operating  the  QSKIP  at  about  9K  should 
provide  very  good  performance  to  about  32p.m  with  a  calculated  NEP  of  about  7x1 0’  Watts. 
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2.3  DETECTOR  FABRICATION 


Fabrication  of  the  hybrid  LTS/HTS  photodetector  with  readout  involved  finding  a 
compromise  between  the  disparate  materials  involved  in  LTS  vs.  HTS  processing.  Possible 
imdesirable  interactions  between  HTS  and  LTS-related  areas  include: 

•  Loss  of  oxygen  from  YBCO  to  other  materials,  such  as  Nb  and  SiOx 

•  Diffusion  of  materials  such  as  Si  or  Nb  into  the  YBCO 

•  Degradation  of  the  Nb  junction  properties  due  to  the  roughness  of  YBCO 

In  practice  the  most  significant  problem  was  deterioration  of  the  YBCO  Tj  due  to  overlayers 
of  SiOx,  which  robbed  oxygen  fi-om  the  YBCO.  Removal  of  the  SiOx  from  the  YBCO  followed  by 
an  anneal  resulted  in  recovery  of  the  T<.  to  within  a  couple  of  degrees  of  the  original  value.  This 
residual  depression  of  Tj.  was  presumably  due  to  some  diffusion  of  Si  into  the  YBCO. 

The  initial  process  for  a  magnetically  coupled  photodetector  and  readout  involved  three  such 
SiOx  depositions,  which  resulted  in  the  T^  being  supressed  to  zero.  Subsequent  designs  for  a  direct- 
coupled  photodetector  and  readout  involved  only  one  such  deposition.  Depending  on  the  particular 
YBCO  film  this  typically  lead  to  supression  of  Tg  by  anywhere  from  zero  to  five  degrees  by  this 
single  step.  In  practice  the  sum  total  of  all  process  steps  resulted  in  suppression  of  T^  by  as  much  as 
twenty  degrees,  a  serious  degree  of  degradation. 

It  should  be  noted  that,  in  retrospect,  an  all-epitaxial  approach,  using  only  epitaxial 
insulators  compatible  with  YBCO,  would  lead  to  less  process-related  damage  to  the  YBCO  than  the 
hybrid  LTS/HTS  approach  pursued  on  this  program.  Needless  to  say,  this  was  not  obvious  at  the 
outset  of  the  program. 
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The  remainder  of  this  section  discusses  the  processes  developed  for  both  magnetically 
coupled  and  direct  coupled  HTS/LTS  photodetector/readout  structures. 

2.3.1  Magnetically  Coupled  Detector 

The  initial  approach  to  the  hybrid  HTS/LTS  photodetector  was  to  have  a  YBCO 
photodetector  on  the  bottom,  then  a  layer  of  SiO^  to  isolate  the  photodetector  from  an  overlayed  Nb 
SQUID  readout  circuit,  which  was  in  turn  overlayed  with  another  SiO^  layer  to  provide  isolation 
from  the  top  Au  layer,  which  served  a  a  light  shield  and  contact  layer.  This  approach  is  outlined  in 
Figure  2.3,  which  also  details  the  individual  process  steps  involved. 


10)  Au  contact,  flux  bias 
&  light  shield 
deposition  &  liftoff 

9)  Insulator  deposition 
&  etch  (Si02) 


Nb  wire  up  deposition 
&  liftoff 

Mo/Ti  resistor  depositio 
&  liftoff 

Si02  deposition  &  RIE 
Nb  JJ  RIE 

Nb/AlOx  trilayer  liftoff 


3)  Insulator  deposition 
&  etch  (Si02) 


Au  contact  deposition, 
liftoff,  &  anneal 
YBCO  deposition 
&  lift  off 

Buffer  layer  deposition 
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Figure  2.3  Schematic  process  sequence  for  HTS  meander  photodetector  with  magnetically  coupled 
LTS  SQUID  readout.  In  practice  the  SQUID  readout  is  accomplished  with  a  series  array  of 
SQUIDs. 
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Two  versions  of  this  photodetector  were  fabricated  on  the  wafer:  one  so-called  “single-slit”, 


I 

which  had  a  simple  loop  of  YBCO  with  a  series  array  of  Nb  SQUIDs  overlayed  on  the  slit,  and  a 
“meander”  version  in  which  the  SQUIDs  were  layed  over  the  central  slit  of  the  meander.  Mask 
drawings  for  the  meander  version  are  shown  Figure  2.4,  while  Figure  2.5  shows  a  detail  of  the 
I  SQUID  array  on  the  single-slit  version. 


Au  light  shield/flux  bias 


Completed  sensor 


Figure  2.4  Excerpts  from  mask  drawings  for  a  YBCO  meander  photodetector  with  series  array  of 
Nb  SQUIDs  for  readout. 
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One  of  the  concerns  at  the  start  of  the  program  was  whether  the  quality  of  the  Nb  junctions 
would  be  compromised  by  the  relatively  rough  surface  of  the  YBCO  below.  This  was  not  however 
a  major  problem. 

From  Figure  2.3  one  can  see  that  a  total  of  three  layers  of  SiO^  were  required  in  this 
process,  including  one  required  for  defining  contact  windows  to  the  Nb  junctions.  Each  deposition 
of  SiOx  degraded  the  YBCO,  resulting  in  it  being  non-superconducting  by  the  time  the  process  was 
completed. 


4-|jm  detector  slit  6-pm  JJ  Mo  resistor  Au  flux-bias 


Figure  2.5  Detail  of  the  mask  layout  for  a  magnetically  coupled  “single  slif’  photodetector  and 
readout,  showing  three  of  nine  SQUlDs  in  series.  The  Au  flux-bias  line  meanders  back  and  forth 
over  the  SQUIDs,  and  is  used  to  set  the  desired  operating  point  of  the  SQUIDs. 
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2.3.2  Direct  Coupled  Detector 


Partly  because  of  the  problems  with  SiO^,  and  because  of  a  general  desire  to  simplify  the 
process,  a  new  photodetector  concept  was  invented  —  the  direct-coupled  photodetector,  where  the 
Nb  SQUID  makes  direct  electrical  contact  to  the  YBCO  through  a  layer  of  in-situ  Au.  The  layer  of 
Au  was  necessary  in  order  to  ensure  that  the  Nb  would  not  steal  oxygen  from  the  YBCO.  Such 
direct  Nb-YBCO  contacts,  without  Au,  have  extremely  high  contact  resistance,  of  order  lO'"*  Q-cm 
or  higher. 

The  fabrication  sequence  for  the  direct-coupled  photodetector  is  illustrated  Figure  2.6.  The 
contact  area  between  the  Nb  SQUID  and  YBCO  photodetector  was  deliberately  made  large  (more 
than  2000  pm^)  in  order  to  obtain  a  high  enough  critical  current  for  the  contact. 


8)  Au  contact  &  light  shield 
deposition  &  lift  off 


7)  Nb  wiring  deposition 
&  lift  off 

|6)  Mo/Ti  resistor  deposition 
&  lift  off 
5)  Si02  deposition 
&  RIE 
|4)  Nb  junction  RIE 

3)  Nb/AlOx  trilayer  liftoff 


2)  Contact  deposition 
&  anneal 

1)  YBCO  deposition 
&  mill 


.a:!; 


Figure  2.6  Schematic  process  sequence  for  HTS  meander  photodetector  with  direct-coupled  LTS 
SQUID  readout.  Only  one  SiO^  layer  is  required  in  this  process,  in  order  to  make  contact  windows 
to  the  Nb  junctions. 
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Figure  2.7  Section  of  the  mask  drawing  for  the  direct-coupled  photodetector,  showing  the  meander 
path  with  light  shield  on  one  branch,  and  light  opening  (through  the  SiO^)  on  the  other  branch,  and 
a  detail  of  the  Nb  SQUID  and  ground  contact  area. 
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3.0  EXPERIMENTAL  APPROACH  AND  MEASUREMENTS 

In  this  section,  we  describe  the  experimental  approach  used  and  the  measurements 
made.  All  the  equipment  and  engineering  effort  used  to  design  the  experimental  test 
station  were  done  on  Westinghouse  funds.  The  experimental  test  station,  shown  in  Figure 
3.1,  has  a  generic  capability  for  making  photoresponse  measurements  and  its  use  is 
intended  for  use  on  several  different  programs  as  a  spectrophotometer  for  measuring 
photoresponse  of  superconducting  and  semiconducting  photophotodetectors.  The 
operation  of  the  test  station  and  date  recording  is  computer  controlled  using  the  LabView 
for  Windows  software  package.  The  optics,  designed  to  operate  over  a  very  vdde  spectral 
band  (3 pm- 100pm),  are  described  in  section  3.1.  The  cryogenics  of  the  test  station  is 
designed  to  operate  between  5-300K  and  is  described  in  section  3.3.  Such  operating 
characteristics  are  compatible  with  testing  superconducting  quantum  photodetectors 
under  development  on  this  program.  Part  of  the  experimental  test  station’s  electrical 
subsystem  was  specifically  designed  for  interface  with  a  SQUID  readout  from  QSKIP, 
and  this  is  described  in  section  3.2.  A  special  packaging  approach  for  handling  the  QSKIP 
was  developed  and  this  is  described  in  section  3.4.  Finally,  in  section  3.5.  we  describe  the 
experimental  results  on  the  QSKIP. 

3.1  OPTICAL  SETUP 

The  optical  setup,  assembled  on  an  aluminum  breadboard  table,  includes:  (1)  a 
monochrometer,  (2)  reflecting  optics,  (3)  cooled  long  pass  filters,  (4)  optical  windows,  (5) 
thermopile  reference  sources,  (6)  black  body  source,  and  (7)  chopper.  A  diagram  of  the 
optical  setup,  including  the  black  body  source,  focusing  optics,  and  monochrometer  is 
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shown  in  Figure  3.2.  Also  included  in  this  diagram  are  the  location  of  the  optical 
elements  in  the  dewar.  The  optical  elements  in  the  dewar  include:  (1)  a  window,  (2) 
cooled  (77K)  filters  (and  attenuators),  (3)  gold  coated  mirrors  for  sampling  the  incident 
and  reflected  optical  beam  onto  thermopiles  (not  shown)  at  77K,  and  (5)  the  sample 
location.  The  sample  is  located  on  a  cold  finger  whose  temperature  can  be  varied  from  5- 
300K.  A  description  of  the  components  used  and  their  operation  follows. 

The  black  body  source  is  made  up  of  a  tungsten  filament  operated  at  2400K.  The 
filament  selected  is  fi'om  an  OSRAM  Halogen  projector  lamp  part  number  HLX  64655 
EHJ.  This  lamp  is  rated  to  operate  at  3400K  and  puts  out  lO'*  lumens  for  an  average  life  of 
50  Hours.  The  filament  is  enclosed  in  a  quartz  envelope  and  operates  fi-om  a  24  DC 
power  supply.  To  use  this  lamp  as  a  black  body  source,  we  first  removed  the  quartz 
envelope  because  it  does  not  transmit  much  beyond  4pm.  The  tungsten  filament  was 
enclosed  in  a  specially  made  vacuum  box  that  had  provisions  for  replaceable  optical 
windows  and  powering  of  the  filament.  The  replaceable  optical  windows  included  KRS-5 
(0.6pm-  40pm),  clear  ZnSe  (0.3pm-  18pm),  and  single  crystal  silicon  (1.1pm-  300pm, 
with  a  blocking  band  between  12  pm-  40pm).  With  these  windows  a  very  broad  spectral 
band  could  be  covered  and  in  particular  the  3  pm-  100pm  spectral  band  we  are  interested 
in. 

To  maintain  a  longer  life,  the  tungsten  filament  was  operated  at  a  lower  voltage 
than  the  rated  24  volts.  The  tungsten  filament  was  enclosed  in  the  vacuum  box  and 
evacuated.  Next,  the  vacuum  box  was  backed  filled  with  an  inert  gas  to  a  pressure  of 
about  30  psi.  This  provided  two  advantages:  (1)  the  window  coating  rate  from  the 
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tungsten  filament  was  reduced,  and  (2)  the  filament  could  be  run  hotter.  Out  of  Argon, 
Krypton  and  Xenon  we  selected  Krypton  because  Xenon  was  too  expensive.  With  this 
design  the  filament  and  the  optical  windows  could  be  readily  replaced. 

Only  Reflective  optics  are  used  to  focus  the  black  body  radiation  onto  the 
monochrometer  entrance  slit.  Two  gold  coated  90°  off-axis  paraboloidal  mirrors  (part  # 
02  POA  013),  fi'om  Melles  Griot,  collimate  and  refocus  the  black  body  radiation  onto  the 
monochrometer’s  entrance  slit.  A  metal  blade  chopper,  located  less  than  3mm  in  fi-ont  of 
the  monochrometer  entrance  slit,  when  rotated  modulates  the  black  body  radiation.  The 
radiation  entering  the  monochrometer  is  filtered  by  a  grating  before  it  emerges  at  the 
monochrometer’s  exit  slit. 

Six  gold  coated  gratings  spectrally  filter  and  provide  radiation  within  the  3  pm- 
100pm  spectral  band  of  interest.  The  six  gratings  Eire  assembled  on  two  turrets  with  each 
turret  holding  three  gratings.  The  monochrometer,  model  SpectraPro-275  from  Acton 
Research  Corporation,  holds  one  three  grating  turret  at  a  time.  Each  turret  is  easily 
installed  in  the  monochrometer  and  no  realignments  are  required.  This  six  grating  set 
covers  the  3pm-  100pm  spectral  band  as  listed  in  Table  3.1.  Turret  number  1  contains 
gratings  (1)  through  (3)  spanning  3-30pm.  Turret  number  2  contains  gratings  (4)  through 
(6)  spanning  30-1 00pm.  The  monochrometer  is  used  to  select  a  narrow  wavelength  band 
(with  higher  order  harmonics)  from  the  spectrally  broad  band  black  body  radiation,  and 
this  optical  signal  is  collimated  at  the  monochrometer  exit  slit  with  a  gold  coated  90°  off- 
axis  paraboloidal  mirrors  (part  #  02  POA  013),  from  Melles  Griot.  The  collimated  optical 
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signal  is  directed  at  the  photodetector  through  an  optical  window  in  the  dewar,  see  Figure 
3.3. 


# 

Grooves  /mm 

Blaze 

wavelength 

Acton 
part  # 

Operating 

Band 

AX  (for  3mm  slit) 

1 

150  g/mm 

4pm 

1-015-4 

3-6pm 

72nm 

2 

75  g/mm 

10pm 

1-075-10 

6-1 6pm 

144nm 

3 

20  g/mm 

22.5pm 

1-02-22.5 

16-30pm 

540nm 

4 

20  g/mm 

45pm 

1-02-45 

30-55pm 

540nm 

5 

13.3  g/mm 

67pm 

1-013.3-67 

55-85pm 

810nm 

6 

7.9  g/mm 

112pm 

1-07.9-112 

85-168pm 

1365nm 

Table  3.1.  The  list  of  gratings  and  their  operating  characteristics  used  in  this  program. 
The  first  three  gratings  were  mounted  on  turret  number  1  and  the  next  three  grating  were 
mounted  on  turret  number  2. 


Part  Number 

Cut-on  X 

Transmission  @  X  Cut-on 

PN  LP-3000-F 

3pm 

5% 

PN  LP-6000-F 

6pm 

5% 

4-8pm  P-12 

12pm 

20% 

5-lOpm  P-16 

16pm 

35% 

8-16pm  P-22 

25pm 

20% 

15-25pm  P-40 

40pm 

20% 

30-40pm  P-62 

62pm 

20% 

Table  3.2.  List  of  long  pass  filters  use  in  the  optical  set-up.  Filters  number  1  and  2  were 
purchased  from  the  Infrared  Optical  Products,  Inc.  Filters  number  3, 4, 5, and  6  were 
purchased  fi'om  Infrared  Laboratories,  Inc. 

Inside  the  dewar,  see  Figure  3.3,  the  collimated  optical  is  filtered  to  remove  the 
higher  order  harmonics.  Seven  long  pass  filters  are  needed  to  operate  over  the  3pm- 
100pm  spectral  band,  and  these  are  listed  in  Table  3.2. 

The  filters,  1”  in  diameter,  are  mounted  in  an  eight-position  wheel  and  cooled  to 
77K.  With  a  knob,  mounted  outside  the  dewar,  we  place  one  of  the  seven  filters  (or  no 
filter)  in  the  optical  beam  path.  Cooling  these  filters  to  77K  greatly  reduces  their  black 
body  radiation.  A  second  wheel,  (cooled  to  77K  and  located  in  front  of  the  filter  wheel 
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and  inside  the  dewar)  contains  neutral  density  attenuators  that  pass  10%  and  1%  of  the 
incident  optical  radiation.  The  second  wheel  is  independently  controlled  with  a  second 
knob  mounted  outside  the  dewar.  This  two  wheel  arrangement  allows  us  to  filter  out  the 
higher  order  harmonics  and  attenuate  the  narrow-band  optical  signal.  The  combination  of 
filters  and  gratings  used  for  spanning  the  3p,m-  100pm  spectral  band  is  given  in  Table 
3.3.  The  last  column  indicates  the  location  of  the  filters  in  the  wheel  so  the  knob  moimted 
outside  the  dewar  can  be  used  to  set  the  dial  at  the  appropriate  number  (or  filter). 
Similarly,  the  description  of  the  wheel  with  the  attenuators  is  given  in  Table  3.4  and  the 
knob’s  dial  settings  are  listed  in  the  last  coliunn. 


Operating 

Band 

#  Grating  &  Part 

Filter 

Knob  Setting 

Open 

n/a 

None 

240 

3-6pm 

(1)  1-015-4 

PN  LP-3000-F 

340 

6-12pni 

(2)  1-075-10 

PN  LP-6000-F 

440 

12-16fj.m 

(2)  1-075-10 

4-8pm  P-12 

540 

16-26pm 

(3)  1-02-22.5 

5-lOpm  P-16 

640 

26-30pm 

(3)  1-02-22.5 

8-16)j.m  P-22 

740 

30-45pn[i 

(4)  1-02-45 

8-16pm  P-22 

740 

45-55pm 

(4)  1-02-45 

15-25pm  P-40 

840 

55-70pm 

(5)  1-013.3-67 

15-25pni  P-40 

840 

70-85pm 

(5)  1-013.3-67 

30-40pm  P-62 

940 

85-lOOpm 

(6)  1-07.9-112 

30-40pin  P-62 

940 

X  XX  WWX  MXXV4.  V^X«XXXXX^U  XV^X  V/|^WX  IXI.XXX^  Xii  X  Vfvyjcxxxx  O^WV/XX  MX  X/UXXVX. 

Once  the  operating  band  is  selected,  the  corresponding  grating  is  rotated  into  place  in  the 
monochrometer  and  the  appropriate  filter  in  the  dewar  is  rotated  into  the  beam  path  with 
the  knob  on  the  dewar. 


After  being  spectrally  filtered  inside  the  dewar,  the  filtered  optical  signal  is 
divided  into  three  parts.  One  part  is  reflected  by  a  gold  mirror  and  is  used  as  a  reference. 
A  second  part  is  reflected  from  a  thin  film  held  at  the  same  temperature  as  the 
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photodetector.  The  thin  film  is  the  same  material  that  the  photodetector  is  made  from. 


The  third  part  of  the  filtered  optical  signal  reaches  the  photodetector  located  on  the  cold 
finger. 

Two  thermopiles  (operated  at  77K  and  mounted  on  a  cross  member  located  inside 
the  dewar’s  radiation  shield,  shown  in  Figure  3.3),  are  used  to  measure  the  intensity  of  the 
incident  and  reflected  beams.  With  the  measured  intensities  of  the  reflected  and  incident, 
the  detector’s  normalized  photoresponse  is  computed.  Both  thermopiles  are  model  2M, 
manufactured  by  Dexter  Research  Center.  Improved  thermopile  responsivity  is  realized 
because  they  are  operated  in  a  vacuum  and  at  77K.  The  responsivity  of  both  thermopiles 
have  been  calibrated  with  a  HeNe  laser.  The  responsivity  of  the  thermopile  used  to 
measure  the  incident  filtered  optical  signal  is  1425.5  Volts/Watt.  The  responsivity  of  the 
thermopile  used  to  measure  the  reflected  optical  signal  from  the  photodetector  is  1 104.19 
Volts/Watt.  The  thermopile  output  signals  are  amplified  and  recorded  on  a  PC.  The 


electrical  setup  is  described  in  the  next  section. 


Operating  Band 

Function 

Characteristics 

Knob  Setting 

3-100|i.m 

Open 

Open 

910 

3-lOOpm 

Attenuator 

1%  Pass 

810 

3-100|im 

Attenuator 

10%  Pass 

710 

3-100|j.m 

Open 

Open 

610 

3-lOOpni 

Open 

Open 

510 

3-1 00pm 

Open 

Open 

410 

3-lOOpm 

Open 

Open 

310 

3-lOOpm 

Open 

Open 

210 

Table  3.4.  Arrangement  of  attenuators  in  the  second  wheel.  Out  of  the  eight  positions  in 
this  wheel  only  two  are  used,  for  the  1%  and  10%  neutral  density  attenuators.  With  the 
two  wheel  arrangement,  these  attenuators  can  be  placed  in  series  with  the  filters  to  vary 
the  intensity  of  the  filtered  optical  beam. 
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Figure  3.1.  Shown  is  the  experimental  test  station  for  measuring  photodetector  photoresponse.  The  dewars  shown  in  the  background 
are  used  to  cool  the  radiation  shield  to  77K  and  the  cold  finger  as  low  as  5K.  Also  visible  are  the  electronics,  dewar,  optics  and 
computer  control. 
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Figure  3.2.  The  layout  of  the  optical  system  is  shown,  including:  the  monochrometer,  the  paraboloidal  mirrors,  and  the  location  of  the 
dewar  components. 


Figure  3.3.  Details  of  the  dewar  design  are  shown,  with  the  extra  liquid  nitrogen  tank  to  cool  the  radiation  shield  to  77K.  The  cold 
finger  is  cooled  by  Helium  to  as  low  as  5K.  The  dual  filter  wheel  assembly  is  attached  to  the  top  of  the  radiation  shield.  A  cross 
member  below  the  dual  filter  wheel  assembly  contains  the  two  thermopiles  used  to  measure  the  incident  and  reflected  optical  beams. 


3.2  ELECTRICAL  SETUP 


The  electrical  setup  included  provisions  for:  (1)  PC  computer  data  acquisition,  (2) 
signal  processing  electronics,  (3)  photodetector  bias  circuits,  (4)  electrical  interface  to 
optical  elements,  and  (5)  temperature  control  and  monitoring.  A  significant  effort  was 
required  for  the  electrical  subsystem  and  this  is  described  next. 

PC  computer  data  acquisition  is  used  to  control  the  entire  experimental  setup.  A 
PC  with  a  486  microprocessor  is  used  to  operate  the  equipment  through  the  GPIB. 
Incorporated  into  the  PC  are  two  cards  made  by  National  Instruments:  (1)  an  AT-GPIB 
interface  for  IBM  PC/AT  EISA  Bus  PC  including  NI-488.2,  and  (2)  a  multifunction 
Analog  and  Digital  I/O  Bus  for  IBM  PC,  model  AT-MIO-16L-25.  The  computerized  data 
acquisition  software  is  based  on  National  Instruments  LabView  for  Windows.  All  the 
electronics  are  initialized  through  the  PC  computer  interface,  with  provisions  for  a 
manual  override.  The  initial  settings  have  been  programmed  into  the  Virtual  Instruments 
(VI)  formed  with  the  LabView  software.  The  VI  incorporates  all  the  electronic 
instruments  on  the  PC  CRT  screen  and  provides  control  knobs  for  operating  these 
through  an  GPIB.  Once  the  photodetector  sample  has  been  cooled  and  properly  biased, 
the  data  recording  is  made  with  the  PC  computer.  The  PC  computer  reads  many  different 
parameters  from  the  instruments  and  these  are  listed  in  Table  3.5.  The  reason  for  all  these 
reading  follow. 

The  first  three  rows  are  used  to  establish  the  AC  incident  photon  flux.  The 
incident  flux  is  measured  with  a  thermopile  when  the  chopper  blade  is  open  and  closed. 
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The  difference  between  these  reading  gives  us  the  change  in  photon  flux  and  this  is  used 


to  compute  the  AC  photocurrent  response.  The  AC  measurement  is  used  to  remove  the 
photon  background  which  is  very  significant  in  the  3- 100pm  spectral  band. 


# 

Label 

Purpose 

1 

High  Incident  Photon  Reading 

Incident  flux  recording  with  chopper  blade  open 

2 

Low  Incident  Photon  Reading 

Incident  flux  recording  with  chopper  blade  closed 

3 

Delta  1-2  COMPUTATION 

Incident  flux  AC  amplitude  computation 

4 

High  Reflected  Photon  Reading 

Reflected  flux  recording  with  chopper  blade  open 

5 

Low  Reflected  Photon  Reading 

Reflected  flux  recording  with  chopper  blade  closed 

6 

Delta  4-5  COMPUTATION 

Reflected  flux  AC  amplitude  computation 

7 

Detector  High  Voltage  Reading 

Voltage  across  detector  with  chopper  blade  open 

8 

Detector  Low  Voltage  Reading 

Voltage  across  detector  with  chopper  blade  closed 

9 

Delta  7-8  COMPUTATION 

Voltage  change  across  SC  should  remain  ZERO 

10 

Detector  High  Bias  Current 

Measures  bias  current  with  chopper  blade  open 

11 

Detector  Low  Bias  Current 

Measures  bias  current  with  chopper  blade  closed 

12 

Delta  10-11  COMPUTATION 

Insures  bias  current  remains  the  same 

13 

Detector  Temperature 

Reads  the  detector’s  temperature 

14 

SQUID  Temperature 

Reads  the  cold  finger’s  temperature 

15 

Cold  Shield  Temperature 

Reads  the  radiation  shield’s  temperature 

16 

Thermopile  Temperature 

Reads  the  thermopile’s  temperature 

17 

SQUID  Bias  Current 

Ic  SQUID  bias  current 

18 

SQUID  High  Current  Reading 

Current  used  to  balanced  SQUID  with  blade  open 

19 

SQUID  LOW  Current  Reading 

Current  used  to  balanced  SQUID  with  blade  closed 

20 

Delta  18-19  COMPUTATION 

AC  photocurrent  amplitude 

21 

Photon  Wavelength 

Provides  spectral  band  of  incident  photon  flux 

22 

DC  Lock-in  Null  Signal 

Error  signal,  zero  if  18  &1 9  are  adjusted  OK 

23 

SQUID  Flux  Bias 

Net  SQUID  bias  including  Lock-in  error  signal 

24 

Lamp  Voltage 

DC  Bias  voltage  on  Black  Body  lamp  source 

Table  3.5.  List  of  readings  and  computations  produced  with  the  computerized  data 
acquisition  system  used  to  measure  the  photodetector’s  photoresponse  and  tabulated  in  24 
rows. 


Rows  4  through  6  are  used  to  establish  the  AC  reflected  photon  flux.  The  reflected 
photon  flux  is  measured  with  a  thermopile  with  the  chopper  blade  open  and  closed.  The 
reflected  photon  flux  is  normalized  by  the  corresponding  incident  photon  flux.  The 
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difference  between  these  readings  should  yield  the  percentage  of  photon  signal  which  is 
reflected.  With  the  AC  reflected  photon  signal  calculations  are  made  to  compute  how 
much  photon  flux  is  absorbed  by  the  photodetector  to  produce  the  observed  photocurrent 
signal. 

Rows  7  through  9  are  used  to  monitor  the  voltage  across  the  superconducting 
photodetector.  Two  readings  are  made,  one  with  the  chopper  blade  open  and  a  second  one 
with  the  chopper  blade  closed.  If  the  superconducting  photodetector  remains  in  the  zero 
resistance  state  (or  current  state)  the  difference  between  these  two  readings  should  be 
zero.  Taking  a  difference  between  two  reading  circumvents  the  issue  of  removing  DC 
offsets  present  in  the  voltage  reading.  If  the  difference  between  the  two  readings  is  not 
zero  it  signifies  that  the  photodetector  is  not  operating  in  the  zero  resistance  state. 

Rows  10  through  12  are  used  to  monitor  the  current  bias  of  the  superconducting 
photodetector.  Two  readings  are  made,  one  with  the  chopper  blade  open  and  a  second  one 
with  the  chopper  blade  closed.  If  the  superconducting  photodetector  remains  in  the  zero 
resistance  state  (or  current  state)  the  photodetector  difference  between  these  two  readings 
should  be  zero.  Taking  a  difference  between  two  readings  circumvents  the  issue  of 
removing  DC  offsets  present  in  the  reading.  The  combination  of  the  voltage  and  current 
readings  is  our  check  to  ensure  that  the  photodetector  remains  in  the  zero  resistance 
superconducting  state  with  different  levels  of  photoillumination. 

Rows  13  through  16  provide  temperature  information  on  the  photodetector  and 
dewar.  Row  13  represents  the  temperature  from  a  diode  thermometer  mounted  inside  a 
leadless  chip  carrier  package  containing  the  photodetector.  Row  14  is  the  temperature  of 
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silicon  diode  sensor  thermally  attached  to  the  cold  finger.  The  temperature  of  the 
radiation  shield  is  recorded  on  line  15  and  we  assume  that  it  corresponds  to  the 
temperature  of  the  thermopiles  which  are  attached  to  the  radiation  shield.  Row  16  has 
provisions  for  adding  a  second  temperature  sensor  for  specifically  measuring  of  the 
thermopiles  temperature,  which  are  bolted  to  the  cross  member  in  the  radiation  shield. 

The  SQUID  bias  current  is  recorded  in  row  17.  The  bias  current  is  manually 
adjusted  for  placing  the  SQUID  in  the  voltage  state  and  at  maximum  sensitivity. 

Rows  18  through  20  monitor  the  needed  offset  current  for  balancing  the 
photodetector  signal  and  placing  the  SQUID  in  the  symmetrical  operating  mode.  The 
offset  current  is  adjusted  with  the  chopper  blade  open  (Row  1 8)  and  the  chopper  blade 
closed  (Row  19).  The  difference  between  these  readings  (Row  20)  represents  the 
photodetector’s  AC  photoresponse  current.  This  photoresponse  is  normalized  with  the 
incident  flux  level  obtained  from  Rows  1-6.  The  number  of  photons  incident  are 
calculated  from  the  wavelength  (Row  21)  and  the  spectrophotometer  spectral  bandwidth 
Ak  given  in  Table  3.1. 

The  total  offset  current  flowing  in  the  SQUID  is  automatically  trimmed  with  a 
Loek-in  feedback  circuit  and  the  trim  offset  current  level  is  given  in  Row  22.  The  total 
offset  current  is  a  sum  of  the  manually  adjusted  offset  current  level  and  the  Lock-in 
contribution,  and  this  is  given  in  Row  23.  The  voltage  used  to  power  the  black  body 
tungsten  lamp  is  listed  in  Row  24. 


PART  III 


49 


It  should  be  evident  that  each  data  point  recorded  with  the  computer  data 
acquisition  system  is  extensive.  Such  a  practice  permits  better  data  analysis  and  proper 
operation  of  the  superconducting  photodetectors. 

The  Signal  Processing  Electronics  consist  of  custom  made  preamplifiers  and  off- 
the-shelf  electronics.  The  custom  made  preamplifiers,  housed  inside  a  breakout  box 
attached  to  the  dewar,  provide  gain  to  improve  noise  immunity  to  outside  interference. 
The  output  from  the  preamplifiers  is  channeled  to  the  computer  data  acquisition  system 
through  off-the-shelf  electronics. 

The  off-the  shelf  signal  processing  electronics  include:  (1)  a  Keithley  182 
nanovolt  meter  with  an  IEEE  488  interface,  (2)  a  Keithley  705  scanner  with  a  7160 
nanovolt  scanner  card  and  IEEE  488  interface,  and  (3)  an  Ithaco  3970  electro  optical 
Lock-in  measuring  system.  These  instruments  receive  the  outputs  from  the  preamplifiers 
connected  to  the  superconducting  photodetector  as  shown  in  Figure  3.4.  The 
photodetector  is  located  inside  the  dashed  lines  labeled  inside  vacuum  jacket  in  Figure 
3.4. 

The  photodetector  bias  circuits  are  illustrated  in  Figure  3.4  and  the  details  are 
shown  in  Figure  3.5.  A  DC  current  I(QSKIP)  biases  the  photodetector.  The  signal 
developed  is  sensed  with  a  SQUID  located  at  the  photodetector  base.  The  SQUID  is 
biased  with  a  DC  current  Ic  (SQUID).  The  construction  of  the  photodetector  and  the 
readout  SQUID  are  symmetrical  to  provide  a  differential  output  between  the  illuminated 
and  covered  portions  of  the  photodetector.  With  such  a  differential  operation,  the  SQUID 
bias  circuits  provide  bi-directional  flux  bias  (see  Figure  3.5)  that  apply  positive  and 
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flux  null  from  labview 
DAC100 


Figure  3.4.  Diagram  of  preamplifier,  Lock-in  and  photodetector  bias  circuits  used  for  making  photoresponse  measurements. 


negative  currents  to  different  ends  of  the  SQUID.  The  combination  of  offsetting  positive 
and  negative  SQUID  bias  currents  minimizes  ground  currents.  The  bi-directional  DC  and 
AC  flux  modulating  currents  are  produced  through  two  10.2KQ  resistors.  The  amplitude 
and  phase  of  the  SQUID  flux  nulling  currents  are  adjusted  to  be  equal  in  amplitude  and 
180°  out  of  phase  with  a  variable  resistor  and  capacitor  (see  Figure  3.5).  The  flux  nulling 
circuits  sum  three  inputs.  These  inputs  are  the  AC  flux  modulation  signal,  the  V(+)  null 
signal  from  the  Lock-in  and  the  manually  controlled  DC  SQUID  flux  bias  adjust  from  the 
Lab-view  DAC.  Since  the  photodetector  is  DC  coupled  to  a  high  gain  preamplifier 
circuit,  provisions  have  been  included  for  manually  removing  any  DC  offsets  that  may 
saturate  the  preamplifier.  An  “offset  adjust”  current  through  a  20Q  resistor  provides  a 
small  voltage  at  the  preamplifier  inputs  to  null  out  any  DC  offset.  The  preamplifiers  are 
directly  coupled  to  the  SQUID  to  facilitate  characterization  of  the  SQUID’s  I-V  curves. 
The  SQUID  preamplifier  is  an  audio  amplifier  with  a  spectral  voltage  noise  density  of 
less  than  InV/VHz. 

The  electrical  interface  to  optical  elements  includes:  the  monochrometer,  the 
chopper,  the  black  body  tungsten  filament,  and  the  thermopiles.  The  Acton  Research 
monochrometer  has  an  electronic  control  module  that  interfaces  through  an  IEEE  488  bus 
with  the  PC  computer  to  record  the  operating  wavelength.  The  HMS  light  beam  chopper 
220  interfaces  to  the  computer  through  the  National  Instruments  AT-GPIB  interface  for 
IBM  PC/AT  EISA  Bus  PC  including  NI-488.2.  The  Black  body  tungsten  filament  is 
powered  by  a  remote  controlled  Sorensen  SRL  20-12  power  supply  (0-25VDC,  0-15 
Amp).  The  Sorensen  power  supply  is  controlled  manually  or  with  an  automatic  option. 
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The  output  of  the  supply  is  adjusted  to  a  maintain  a  fixed  photon  flux  level  as  determined 
by  the  thermopiles  inside  the  dewar.  The  thermopiles  are  connected  to  special 
preamplifier  circuits  with  a  gain  of  11  (see  Figure  3.6)  and  are  calibrated  to  read  the 
incident  photon  flux  power.  To  extend  the  thermopiles’  frequency  response,  the  coaxial 
shield’s  effective  shimting  capacitance  is  reduced  by  the  preamplifier  circuits. 

The  temperature  control  and  monitoring  is  centered  around  two  Lake  Shore 
93CA  temperature  controllers  connected  via  an  IEEE  488  bus  to  a  PC.  Each  controller 
accepts  inputs  from  two  temperature  sensors  and  can  power  one  heater.  The  heater  output 
in  combination  with  either  temperature  sensors  provides  for  temperature  control.  The 
heater  output  is  a  proportional  analog  signal  (not  digital)  to  minimize  heater  induced 
electrical  interference.  The  temperature  sensors  are  calibrated  silicon  diodes  with  a 
temperature  range  between  2K  and  300K.  The  sensor’s  calibration  curves  are  entered 
inside  the  93CA  to  provide,  at  low  temperatures,  a  resolution  better  than  0.  IK. 

3.3  CRYOGENIC  SETUP 

The  dewar  used  in  this  setup  was  a  model  MTD  150  made  by  Lake  Shore 
Cryotronics  (see  Figure  3.3).  The  dewar  includes  100  coaxial  electrical  inputs  for 
electrical  access  to  the  packaged  sample.  The  dewar  includes  a  socket  designed  to  accept 
packaged  devices  in  a  leadless  chip  carrier  (LCC)  84  pin  package.  A  square  hole  is 
machined  inside  the  socket  to  accommodate  an  aluminum  pedestal  which  is  attached  to 
the  cold  finger.  Several  modification  were  made  to  this  dewar  and  they  are  described. 

The  original  dewar  arrangement  made  use  of  a  single  clamp  to  make  electrical  and 
thermal  contact  between  the  LCC  package  and  the  dewar’ s  socket  and  aluminum 
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SQUID  flux  bias 


of:  (1)  bias  circuits  used  for  biasing  the  photodetector  and  the  SQUID,  (2)  the  SQUID  preamplifier  readout  circuits, 
monitoring  the  SQUID  critical  current  I^  and  the  flux  bias. 
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pedestal.  We  have  incorporated  two  independent  pressure  clamps  to  electrically  and 
thermally  contact,  respectively,  the  LCC  package  to  the  socket  and  the  aluminum 
pedestal.  The  LCC  socket  is  soldered  to  a  small  printed  circuit  board  which  is  attached  to 
a  second  fan-out  board  with  flexible  thin  wires.  The  small  printed  circuit  board  has 
provisions  for  attaching  a  clamp  that  clamps  down  the  LCC  package  inside  the  socket. 
This  clamp  is  independent  of  a  second  clamp  which  is  used  to  press  the  package/socket 
assembly  onto  the  aluminum  pedestal  and  make  good  thermal  contact.  With  two  spring 
loaded  clamps  good  electrical  and  thermal  contact  are  independently  established.  This  is 
important  given  the  fact  that  the  packaged  devices  are  loaded  into  the  dewar  at  room 
temperature  and  cycled  between  5-300K.  Cycling  devices  over  such  a  wide  temperature 
range  requires  provisions  for  accommodating  mechanical  mismatches  in  the  thermal 
expansion  coefficients.  The  two  clamps  provide  the  additional  degree  of  freedom  to 
accomplish  this  feature. 

The  modified  dewar  included  two  separate  cryogenic  inputs.  The  radiation  shield, 
the  two  filter  wheels,  and  the  thermopiles  were  cooled  by  liquid  nitrogen  through  separate 
cryogenic  lines.  This  permitted  us  to  cool  most  of  the  dewar  mass  with  liquid  nitrogen 
and  saved  both  cool  down  time  and  liquid  helium  consumption.  The  second  cryogenic 
input  was  fi'om  liquid  helium  that  allowed  us  to  cool  the  sample  to  about  5K.  In 
operation,  liquid  nitrogen  and  liquid  helium  were  transferred  at  the  same  time  to  speed  up 
cool  down.  After  cool  down,  the  liquid  nitrogen  valve  was  turned  off  and  the  radiation 
shield  remained  close  to  77K.  The  cold  finger  temperature  was  varied  between  5-77K. 
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To  reduce  electrical  interference,  a  Mu-metal  magnetic  shield  was  incorporated  in 
the  dewar.  The  shield,  purchased  from  Advance  Magnetics,  was  custom  designed  and 
made  from  RMA  1084  material.  The  shield  was  assembled  from  several  pieces  so  it  could 
be  incorporated  into  the  dewar.  As  designed  the  shield  was  easily  added  to  the  dewar  after 
all  the  electrical,  cryogenic,  vacuum,  and  optical  elements  were  incorporated. 

As  configured,  the  dewar  was  routinely  used  to  cool  down  and  measure 
photodetectors  down  to  5K.  However,  operation  at  5K  required  reducing  the  thermal  load 
from  the  100  coaxial  lines.  For  testing  single  photodetectors,  25  coaxial  lines  were 
sufficient  and  75  coaxial  lines  were  discormected  from  the  cold  finger  stage  to 
significantly  reduce  the  thermal  load.  The  75  lines  can  be  readily  reattached  by  removing 
three  plastic  inserts. 

3.4.  DETECTOR  PACKAGING 

Testing  of  LWIR  QSKIP  IC  requires  provisions  for  dealing  with  large  thermal 
expansion  mismatch  between  the  photodetector  and  package  material.  Additionally, 
operation  at  liquid  helium  temperatures  requires  excellent  thermal  contact  between  the 
package  and  the  cold  finger.  We  developed  a  method,  on  Westinghouse  fimds,  for 
packaging  superconducting  IC’s  in  modified  ceramic  packages  which  resolve  the  thermal 
expansion  mismatch  and  thermal  drops  across  the  package. 

The  solution  we  have  developed  solves:  (1)  the  thermal  expansion  mismatch,  and  (2) 
thermal  drops  across  the  package.  The  packaging  approach  is  compatible  with 
conventional  bonding  of  integrated  circuits  into  packages  that  fit  in  conventional  sockets. 
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The  socket  used  is  for  an  84  pin  LCC  package,  Kyocera  part  number  PB-0289,  and  is 
integrated  into  the  Lake  Shore  Cryotronics  MTD-150  open  cycle  dewar. 

The  thermal  mismatch  is  resolved  by  using  a  metallic  base  material  (for  good 
thermal  conductivity)  that  has  an  thermal  expansion  coefficient  close  to  the  substrate  the 
superconductor  material  is  grown  on.  YBCO  is  typically  grown  on  LaA103  and  NdGa03 
substrates.  Table  3.6  lists  thermal  expansion  coefficient  for  several  materials.  For  a 
LaA103  substrate,  beryllium,  niobium  and  tantalum  offer  very  good  thermal  matches 
(1300  vs  1300,  1290,  and  1279).  Niobium  and  tantalum  are  potentially  a  good  choices. 
However,  as  niobium  becomes  superconducting  at  9.2K  (vs  4.43K  for  tantalum)  its 
thermal  conductivity  also  begins  to  decrease  drastically  at  9.2K.  For  operating  down  to 
about  4.5K  it  is  better  to  use  tantalum  and  this  is  out  choice  material. 


MATERIAL 

AL/L  X  10* 

MATERIAL 

AL/L  X  10 

LaA103 

1300 

SS  Monel 

2450 

Niobimn 

1290 

SS  301 

2300 

Tantalum 

1270 

SS  321 

2600 

Molybdenum 

894 

SS316 

2770 

Iron  Fe 

1890 

SS304 

2790 

Beryllium 

1300 

Titanium 

1700 

Aluminum 

3910 

Poly-Alumina 

791 

Table  3.6.  Total  thermal  contraction  difference  AL/L  between  300K  to  77K  for  different 
candidate  substrate. 


Incorporating  a  tantalum  base  into  the  84  pin  LCC  alumina  package  requires 
package  modification.  A  1  cm  diameter  hole  was  drilled  in  the  84  pin  LCC  package  to 
remove  the  alumina  that  the  IC  in  normally  mounted  on  (see  3.7).  The  hole  was 
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sufficiently  small  to  be  confined  within  the  package  well  and  thus  not  damage  the 
electrical  lines  inside  the  package.  The  hole  removes  the  bottom  of  the  IC  well  which  is 
replaced  by  the  tantalum  base(see  Figure  3.7). 

1  21 


1  21 

Figure  3.7.  Schematic  of  84  pin  leadless  chip  carrier  package  with  hole  drilled  in  IC 
cavity.  The  hole,  indicated  by  the  circle,  is  1  cm  in  diameter,  which  is  sufficiently  large 
for  the  IC  to  fit  inside. 

The  base  is  made  from  a  thin  (about  0.05cm  thick)  tantalum  sheet  cut  into  a  2.4 
cm  square.  The  2.4  cm  square  size  is  sufficiently  small  to  fit  within  the  electrical  contact 
fi-ame  located  on  the  back  of  the  ceramic  package,  about  2.86  cm  square.  The  processing 
and  assembly  needed  to  incorporate  this  tantalum  base  are  described  next. 

A  tantalum  sheet  is  first  plated  with  about  lOpm  of  copper  on  both  sides  and  then 
cut  into  2.4  cm  small  squares.  The  plated  tantalum  squares  receive  about  120pm  of 
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Indium  on  the  front  side,  for  a  total  base  thickness  of  about  620^m.  The  LWIR  QSKIP 
devices  are  soldering  to  the  indium-coated  tantalum  base.  The  soldering  is  performed  by 
placing  the  LWIR  QSKIP  device  (with  a  gold  coated  back)  on  a  SnInCd  solder  preform 
in  the  center  of  the  indium-coated  tantalum  base.  This  assembly  is  placed  into  an  oven 
and  heated  to  130°C  and  cooled  to  room  temperature.  This  thermal  cycle  solders  the 
LWIR  QSKIP  device  to  the  indium-coated  tantalum  base  with  the  SnInCd  solder  whose 
melting  temperature  is  95°C.  The  tantalum  base  with  the  LWIR  QSKIP  device  is  glued  to 
the  back  of  the  ceramic  package  as  shown  in  Figure  3.8. 


LEADLESS  CfflP  CARRIER 
PACKAGE  WITH  HOLE 
DRILLED  IN  IC  WELL 


TANTALUM  BASE 

Figure  3.8.  Illustration  of  assembly  of  indium-coated  tantalum  base,  with  IC  attached, 
before  being  glued  to  back  of  ceramic  package. 

With  this  assembly,  the  soft  indium,  located  between  the  tantalum  base  and  the 

alumina  LCC  package,  provides  mechanical  compliance  to  accommodate  thermal 
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expansion  mismatch.  The  fast  setting  glue,  (Super  Bonder  495  made  by  Loctite)  is  used 
for  ease  of  assembly.  When  the  LWIR  QSKIP  device,  the  tantalum  base,  and  package  are 
assembled,  ultrasonic  bonding  is  used  to  bond  the  LWIR  QSKIP  to  the  LCC  package 
frame.  The  completed  assembly  is  shown  in  Figure  3.9,  with  the  electrical  contact  pads 
exposed  on  the  back.  With  this  arrangement,  the  electrical  contact  between  the  socket  and 
LCC  package  is  through  the  original  contact  pads,  located  on  the  package’s  back  side. 

SUPERCONDUCTING  IC  MOUNTED 
HOLE  IN  LEADLESS  CHIP  CARRIER 
AND  SOLDERED  WITH  SnInCd  (95C)  SOLDER 
INDIUM  COATED  TANTALUM 


Figure  3.9.  Assembly  of  Tantalum  base,  with  IC  soldered  to  it,  glued  to  the  back  of  a 
ceramic  package  pre-drilled  with  a  hole  in  the  IC  well. 

The  packaging  scheme  has  been  implemented  and  LWIR  QSKIP  devices  have  been 
placed  in  a  dewar  and  cooled  down  to  about  5K.  Unlike  previous  approaches,  the  LWIR 
QSKIP  device  did  not  sheer  off  the  tantalum  base  after  several  cycles  between  300K  and 
5K.  This  is  a  significant  improvement  over  previous  attempts  where  the  QSKIP  was 
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placed  directly  inside  a  84  pin  ceramic  package.  This  technique  is  appropriate  for 
packaging  different  IC  made  with  high  and/or  low  temperature  superconductors. 


3.5  DETECTOR  TESTING  AND  RESULTS 

The  LWIR  QSKIP  fabricated  and  tested  in  this  program  is  represented  by  the 
diagram  shown  in  Figure  2.2  and  the  photograph  shown  Figure  3.10.  In  this  section  we 
start  with  an  analysis  specific  to  this  photodetector  structure  and  follow  with 
experimental  measurements  taken  on  photodetectors  pictured  in  Figure  3.10.  The  LWIR 
QSKIP  did  not  exhibit  photoresponse  to  light  and  we  examine  potential  reasons  for  this. 

The  device  pictured  in  Figure  3.10  includes  a  YBCO  photodetector  with  a 
niobium  readout  SQUID  located  at  the  bottom.  The  photodetector  is  symmetric  with  each 
side  configured  into  a  serpentine  pattern  to  increase  the  kinetic  inductance.  One  side  of 
the  photodetector  is  shielded  with  a  gold  light  shield  and  the  other  side  exposed  to  light. 
The  SQUID  readout  consists  of  two  resistively  (by  IQ)  shunted  Nb/Al203/Nb  Josephson 
tunnel  jimctions. 

In  operation,  the  serpentine  photodetector  is  biased  by  a  DC  current.  The  0.1pm 
thick  YBCO  photodetector  occupies  a  100pm  x  100pm  square.  Each  serpentine  branch  is 
about  700pm  long  and  11pm  wide.  For  this  geometry,  the  resistivity  of  undegraded 
YBCO  lines  should  be  about  300pQ-cm  at  300K.  However,  we  have  measured  more 
than  six  devices  and  each  one  exhibits  a  resistivity  of  1800pQ-cm  at  300K.  This 
resistivity  is  six  times  higher  than  the  expected  nominal  value  of  300pQ-cm.  Also,  with 
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processing  the  YBCO  films’  transition  temperature  degrades  from  90K  to  about  60K  (see 
Figure  3.11).  These  changes  in  resistivity  and  transition  temperature  are  evidence  that  the 
YBCO  LWIR  QSKIP  processing  causes  damage  to  the  YBCO  superconductor. 

We  attribute  the  processing  damage  to  the  YBCO  films  primarily  to  the  silicon 
dioxide  layer  covering  the  YBCO.  The  sputtered  silicon  dioxide  provides  for  electrical 
isolation  between  the  gold  interconnect  lines  and  the  YBCO  film.  The  first  LWIR  QSKIP 
devices  fabricated  had  YBCO  that  was  not  superconducting  down  to  5K.  It  was 
determined  that  the  sputtered  silicon  dioxide  films  were  silicon  rich.  The  excess  silicon  in 
silicon  rich  silicon  dioxide  is  expected  to  leach  oxygen  from  the  YBCO  films. 
Experiment  were  run  to  add  oxygen  to  the  Argon  sputtering  gas  with  the  expectation  of 
realizing  stochiometric  silicon  dioxide.  With  this  approach,  the  processed  YBCO 
remained  superconducting  after  processing  albeit  they  exhibited  a  reduced  transition 
temperature.  We  stripped  the  oxide  from  these  YBCO  films  and  attempted  to  reoxygenate 
these  film.  Thermal  cycles  in  oxygen  ambient  at  400C  did  not  increase  the  YBCO  films 
transition  temperature  to  90K.  The  reduced  transition  temperature  remained  leading  us  to 
conclude  that  some  silicon  has  diffused  into  the  YBCO  sample  and  was  irreversible. 
Careful  tuning  of  the  silicon  dioxide  deposition  improved  the  YBCO  films.  The  latest 
YBCO  films  remained  superconducting  after  all  the  processing,  which  incorporate  the 
niobium  SQUID.  However  with  the  reduced  transition  temperature  and  increased 
resistivity  we  believe  that  the  degraded  superconducting  characteristics  have  also 
produced  a  significant  degradation  in  the  quasiparticle  lifetime  and  hence  no 
photoresponse  was  observed.  A  solution  to  the  YBCO  degradation  is  to  use  epitaxial 
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SrTi03  as  the  insulator  and  passivation  layer  over  the  YBCO  films.  Such  a  process  is 
being  developed  in  the  laboratory. 

The  SQUID  built  into  the  LWIR  QSBCIP  was  from  niobium.  This  was  selected  to 
reduce  the  risk  of  using  a  YBCO  SQUID.  With  the  niobium  SQUID  noise  and 
performance  issues  are  well  understood,  thus  leaving  only  the  YBCO  photodetector  as  an 
unknovm.  The  niobium  SQUID  incorporated  with  the  YBCO  photodetector  operated.  At 
5K,  the  critical  SQUID  junction  current  was  about  0.2mA.  The  SQUID  I-V  curves  are 
shown  in  an  oscillograph  in  Figure  3.12.  The  SQUID  was  modulated  by  a  magnetic  flux 
and  the  modulation  is  shown  in  Figure  3.13.  The  inductance  of  the  SQUID  loop  was 
about  lOpH,  consistent  with  our  design  goals. 

In  conclusion,  the  calculation  reveal  that  the  LWIR  QSKIP  device  should  provide 
excellent  performance  vvith  an  NEP  of  about  10"  Watts  and  operate  over  a  very  wide 
spectral  bandwidth.  For  YBCO,  the  expected  photoresponse  cutoff  wavelength  is 
expected  to  extend  to  about  30-40|am.  On  this  program,  we  have  developed  performance 
models  and  photodetector  structures  which  should  provide  the  base  for  future 
demonstration  of  superconducting  quantum  photodetectors.  With  the  development  of 
improved  fabrication,  and  QSKIP  designs  this  photodetector  will  be  demonstrated.  The 
main  improvement  for  QSKIP  design  is  a  reduction  of  the  magnetic  field  inductance 
relative  to  the  kinetic  inductance.  Both  process  and  design  improvements  are  being 
pursued  and  we  expect  positive  results  in  the  near  future. 
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Figure  3.10.  A  photograph  of  the  LWIR  QSKIP  designed  and  fabricated  in  this  program. 
The  photodetector  is  made  of  symmetric  YBCO  serpentine  left  and  right  portions  with  a 
niobium  SQUID  readout  located  at  the  bottom. 
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Figure  3.11.  Resistance  of  YBCO  photodetector  versus  temperature  after  all  processing 
was  completed.  There  is  a  small  resistance  offset  which  has  shifted  the  R  vs.  T  curve 
upward.  Processing  of  the  YBCO  films  caused  the  superconducting  transition 
temperature  to  decrease  from  about  90K  to  about  60K. 


Figure  3.12.  The  I-V  curve  of  a  niobium  SQUID  built  into  the  LWIR-QSKIP  device, 
measured  at  4.2K.  The  Y-axis  is  200pV/div  and  the  X-axis  is  200pA/div. 
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Figure  3.13.  The  SQUID  V-O  modulation  curves,  with  magnetic  flux  provided  by  an 
injected  control  current.  The  Y  scale  is  50pV/div  and  the  X-scale  is  200pA/div.The 
measured  period  corresponds  to  a  loop  inductance  of  lOpH,  which  was  the  design  value. 
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Abstract.  The  discovery  of  high-temperature  superconductors  (HTS) 
spawned  many  potential  applications,  including  optical  detectors.  Real¬ 
izing  viable  superconducting  detectors  requires  achieving  performance 
superior  to  competing  and  more  mature  semiconductor  detector  tech¬ 
nologies,  and  quantum  detector  technologies  in  particular.  We  review 
why  quantum  detectors  are  inherently  more  sensitive  than  thermal  or 
bolometric  detectors.  This  sensitivity  advantage  suggests  that  for  oper¬ 
ation  at  cryogenic  temperatures,  we  should  be  developing  only  quantum 
superconducting  detectors.  Accordingly,  we  introduce  and  describe  the 
structure  and  the  operation  of  a  superconducting  quantum  detector  with 
a  superconducting  quantum  interference  device  (SQUID)  readout  circuit. 
The  superconducting  quantum  detector,  consisting  of  a  superconducting 
loop,  produces  a  photosignal  in  response  to  photoinduced  changes  in 
the  superconducting  condensate’s  kinetic  inductance.  The  supercon¬ 
ducting  quantum  detector  is  designed  to  operate  only  in  the  supercon¬ 
ducting  state  and  not  in  the  resistive  or  transition  states. 

Subject  terms:  infrared  technology:  superconducting  sensors;  bolometric  detec¬ 
tors:  quantum  detectors. 

Optical  Engineering  33(3),  697-703  (March  1994). 


1  Introduction 

The  discovery  of  high-temperature  superconductivity  (HTS) 
by  Bednorz  and  Muller*  in  1986  evoked  much  excitement 
and  a  rush  to  find  applications  for  the  new  discovery.  One 
of  the  first  proposed  applications  was  for  making  optical 
detectors,  and  infrared  detectors  in  particular.  This  proposed 
application  was  natural  because  infrared  detectors  operate  at 
liquid  nitrogen  temperatures,  just  like  high-temperature  su¬ 
perconductors.  Enomoto  and  Murakami^  made  the  earliest 
photoresponse  measurements  using  granular  BaPbo.7Bio  3O3 
and  reported  encouraging  results.  Soon,  many  more  papers^"^ 
were  published  reporting  on  the  photoresponse  of  different 
detector  structures  made  in  HTS  materials.  However,  the 
photoresponse  mechanism  was  not  clearly  identified  nor  was 
the  ultimate  expected  detector  performance.  After  much  ex¬ 
amination,  it  was  shown^  that  the  reported  photoresponse 
signals  were  consistent  with  a  bolometric  response.  This  was 
a  very  discouraging  turn  of  events  because  the  sensitivity  of 
bolometric  detectors  is  far  less  than  the  sensitivity  of  quantum 
detectors. 

However,  because  of  inherent  limitations  in  semiconductor- 
based  quantum  detectors,  the  search  for  a  superior  super- 
conductor-based  quantum  detector  persists.  Specifically,  the 
main  improvements  over  semiconductor-based  detectors 
sought  are  a  broadband  photoresponse  and  operation  at  very 
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high  speeds.  Typically,  semiconductors  consume  more  power 
at  higher  operating  speed,  thus  the  desire  to  include  high¬ 
speed  signal  processors  with  the  detectors  is  incompatible 
with  cooling  constraints.  Furthermore,  broadband  sensing,  at 
very  long  infrared  wavelengths,  with  semiconducting  quan¬ 
tum  detectors  demands  much  lower  operating  temperatures 
then  77  K.  HTS  offers  two  attractive  potentials:  first,  super¬ 
conducting  electronics  is  the  only  technology  offering  high¬ 
speed  performance  at  low  power,  and  second,  the  small  en¬ 
ergy  gap  in  HTS  materials  offers  the  possibility  of  broadband 
photodetection.  For  these  reasons  the  development  of  HTS 
broadband  superconducting  detectors  is  being  actively 
pursued. 

In  this  paper,  we  report  on  the  progress  attained  toward 
realizing  HTS  superconducting  quantum  detectors.  We  begin 
by  reviewing  the  inherent  sensitivity  advantage  of  quantum 
detectors  over  bolometric  detectors  in  Sec.  1.1.  This  is  fol¬ 
lowed  by  a  description  of  the  limitations  of  the  conventional 
superconducting  detectors,  given  in  Sec.  1.2.  Finally,  Sec.  2 
describes  our  novel  approach  to  realizing  a  superconducting 
quantum  detector. 

1 .1  Sensitivity  of  Boiometers  and  Quantum 
Detectors 

Bolometers  operate  as  thermal  equilibrium  devices,  whereas 
quantum  detectors  operate  as  nonequilibrium  devices.  The 
performance  of  each  detector  category  is  ultimately  limited 
by  the  detector’s  noise  level.  Photon  shot  noise  is  the  ultimate 
sensitivity  limit  of  an  ideal  detector.  For  very  low  photon 
flux  levels,  the  detector’s  noise  level  sets  the  sensitivity  limit. 
Here  we  review  the  operation  of  each  detector  category  and 
provide  expressions  for  the  sensitivity  limit  of  each.  It  be- 


OPTICAL  ENGINEERING  /  March  1994  /  Vol.  33  No.  3  /  697 


BLUZER  and  FORRESTER 


comes  evident  that  quantum  detectors  are  more  sensitive  than 
bolometric  detectors. 

In  bolometers,  photosignals  are  produced  when  the  pho- 
toabsorbed  photons’  energy  is  quickly  transferred  into  ther¬ 
mal  ized  phonons,  thereby  changing  the  lattice’s  temperature. 
The  transfer  process  of  the  photons’  energy  into  the  lattice 
is  very  fast  relative  to  the  bolometric  detector’s  response  time. 
Although  photoabsorption  in  a  bolometer  occurs  via  inter¬ 
mediate  quantum  excitations  of  electrons  or  phonons,  the 
lifetime  of  these  intermediate  excitations  is  very  short  and 
the  excitation  energy  is  quickly  thermal  ized  to  produce  a 
change  in  the  lattice’s  temperature.  Thus,  the  signal  levels  in 
a  bolometric  detector  correspond  to  the  change  in  the  value 
of  the  lattice’s  temperature  T.  Similarly,  the  sensitivity^’®  of 
a  bolometer  is  limited  by  the  thermal  fluctuations  in  the  lat¬ 
tice’s  temperature  AT^.  For  an  object  with  heat  capacity  C 
connected  to  its  surroundings  by  a  thermal  conductance  G, 
the  object’s  temperature  change  AT  is  related  to  its  power 
exchange  W{t)  with  the  surroundings  by 

dAT 

C—+G^T--W{t)  .  (1) 

dr 

Using  this  expression  and  integrating  over  all  frequencies,^ 
the  mean-square  value  of  the  power  flow  from  the  body 
A  can  be  expressed  as 

^W\=AkBT^G^f ,  (2) 

where  A/  is  the  frequency  band  of  interest.  The  smallest  value 
for  the  object’s  power  fluctuations  with  its  surroundings  oc¬ 
curs  for  the  smallest  value  of  conductance  G,  or  the  con¬ 
ductance  of  free  space  where  a  is  the  Stefan- 

Boltzmann  constant  and  Ap  is  the  surface  area  of  the  object. 
Incorporating  both  equations,  the  mean-square  value^’^  of  the 
power  flow  from  the  body  into  free  space  is 

^W\=  \6Apks(jT^^f  .  (3) 

The  resolution  or  sensitivity  of  a  bolometric  detector  is  lim¬ 
ited  to  the  value  of  the  power  fluctuations  with  the  surround¬ 
ings,  and  this  limit  cannot  be  reduced  by  readout  circuits. 

In  quantum  detectors,  photoabsorption  efficiently  chan¬ 
nels  the  photons’  energy  directly  into  the  production  of  quan¬ 
tum  excitations,  which  may  have  relatively  long  lifetimes, 
creating  a  nonequilibrium.  The  nonequilibrium  consists  of 
an  excess  population  of  excited  states  over  the  thermal  equi¬ 
librium  population.  The  signal  and  the  noise  in  quantum  de¬ 
tectors  correspond  to  the  excess  nonequilibrium  population 
and  the  fluctuations  therein,  respectively.  Thus  for  a  given 
lattice  temperature,  different  nonequilibrium  levels  exist,  un¬ 
like  in  bolometric  detectors.  The  sensitivity  of  a  quantum 
detector  is  limited  by  the  larger  of  ( 1 )  the  minimum  electronic 
noise  in  the  detector  AITq£  or  (2)  the  photon  radiation  shot 
noise  AW^^.  These  noise  quantities  (AW^^,  AWqp)  are  ex¬ 
pressed  as  an  average  power  fluctuation  in  the  quantum  de- 
tector^’^  as 


^WQp  =  AkJAf  , 


AWqp  = 


2'TrcX  ^  dX 
Q\^{hd\kpT)  —  1 


(4a) 

(4b) 


where  c  is  the  speed  of  light,  X  is  the  wavelength,  X^  is  the 
cutoff  response  of  the  detector,  and  h  is  Plank’s  constant.  It 
should  be  evident  from  the  integral  in  Eq.  (4b)  that  the  photon 
radiation  shot  noise  in  the  quantum  detector  is  band  limited 
by  the  cutoff  wavelength  X^..  In  bolometric  detectors,  there 
is  no  spectral  band  limiting  and  hence  the  photon  noise  is 
significantly  higher.  Evaluating  at  77  K,  the  noise  power 
expressions  for  bolometric  and  quantum  detectors  (operating 
at  a  30-Hz  frame  rate),  we  find  that  the  noise  power 
in  a  bolometric  detector  is  significantly  higher  than 
the  noise  power  in  a  quantum  detector,  each  detector  being 
10“^  cm^  in  area.  The  bolometer’s  noise  power  (AW^)'^^  is 
calculated  to  be  10  “  W,  for  a  scene  and  a  detector  operating 
temperature  of  77  K.  The  Johnson  noise  power  for  a 

quantum  10“"^  cm^  detector  operating  at  77  K  is  1.4  X  10“  W, 
Additionally,  the  radiation  shot  noise  for  3  to  5  (8  to  12 

|xm)  operation  is  10“*^  W  (10“‘^  W).  Combining  the  ra¬ 
diation  shot  noise  with  the  detector’s  Johnson  noise  in  quad¬ 
rature,  it  becomes  evident  that  the  noise  in  quantum  detectors 
is  much  less  than  the  noise  in  bolometric  detectors.  These 
calculations  clearly  illustrate  that  at  liquid  nitrogen  temper¬ 
atures,  quantum  detectors  are  inherently  more  sensitive  than 
thermal  detectors.  This  preference  for  bolometric  versus 
quantum  detectors  is  independent  of  the  specific  detector 
signal  generation  mechanism,  e.g.,  polarizibility,  resistivity, 
reflectivity,  etc.  Note  that  in  addition  to  minimizing  the  noise 
levels  in  a  detector,  the  signal  levels  need  to  be  maximized 
to  facilitate  the  performance  of  the  readout  circuits.  The  pho¬ 
tosignal  levels  in  bolometric  detectors  are  maximized  with 
minimum  values  for  C  and  G,  which  ultimately  puts  severe 
restrictions  of  the  operating  frame  rate  of  bolometric  detec¬ 
tors.  In  quantum  detectors,  the  photosignal  levels  are  max¬ 
imized  with  maximum  values  for  the  excitation  lifetimes  and 
much  higher  operating  frame  rates  are  possible. 

1 .2  Conventional  Superconducting  Detectors 

Like  semiconductors,  superconducting  detectors  can  be 
grouped  into  two  main  categories:  quantum  and  bolometric 
(or  thermal)  detectors.  Several  different  detectors  have  been 
realized  in  low-  and  high-temperature  superconductors  and 
they  can  be  classified  into  the  categories  mentioned.  Other 
exotic  mechanisms  of  detection,  such  as  photon-induced 
vortex-antivortex  nucleation  and  photon-enhanced  flux-flow, 
have  been  proposed^  but  are  not  well  understood  vis-a-vis 
the  four  detector  categories  enumerated  below. 

1.  Thermal  detectors,  where  the  energy  from  incident 
photons  changes  some  temperature-dependent  prop¬ 
erty  of  the  material,  such  as  resistance.  Transition  edge 
bolometers  have  been  realized  in  low-temperature  su¬ 
perconductor  (LTS)  and  HTS  materials,  and  their  sen¬ 
sitivity  is  limited  by  the  thermal  energy  fluctuations  of 
the  lattice  and  the  photon  flux.  These  detectors  provide 
attractive  performance  especially  for  long-wave 
(>20  jxm)  detection  where  semiconductors  require  liq¬ 
uid  helium  operating  temperatures. 

2.  Josephson-effect-based  detectors,  which  make  use  of 
the  nonlinear  interaction  between  the  incident  photon’s 
electromagnetic  field  and  the  high-frequency  oscillat¬ 
ing  currents  and  voltages  in  a  biased  Josephson  junction 
to  perform  millimeter-wave  mixing  and  direct  detec- 
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tion.  These  devices  are  particularly  useful  for  milli¬ 
meter-  and  submillimeter-wave  detection. 

3.  Quasipanicle  tunneling  devices,  such  as  superconduc¬ 
tor-insulator-superconductor  (S-l-S)  mixers  and  direct 
detectors,  where  the  incident  photons  are  absorbed  by 
quasiparticles,  allowing  tunneling  events  that  would 
otherwise  be  energy  forbidden  to  affect  the  junctions 
current-voltage  (I-V)  characteristics. 

4.  Pair-breaking  detectors,  where  incident  photons  with 
energy  larger  than  the  energy  gap  break  Cooper  pairs, 
leading  to  changes  in  such  macroscopic  parameters  as 
the  superconducting  gap  and  the  kinetic  inductance. 
Because  of  the  difficulty  of  making  S-I-S  structures  in 
HTS,  this  fourth  detector  category  has  received  the 
most  attention  and  is  described  in  the  following. 


The  fourth  category  is  the  least  developed  from  a  historical 
and  practical  standpoint.  It  can  provide  thermal  as  well  as 
quantum  detectors  critically  depending  on  the  material  pa¬ 
rameters  such  as  quasiparticle  recombination  time,  which  is 
not  relevant  to  the  other  categories  described.  In  HTS  ma¬ 
terials,  the  quasiparticle  lifetime  and  the  theory  for  predicting 
it  are  unknown.  What  is  accepted  is  the  existence  of  Cooper 
pairs  in  HTS  with  a  minimum  energy  2A  required  to  break 
them.  Thus,  a  basic  mechanism  of  a  superconducting  quan¬ 
tum  detector  is  to  have  photons  with  energy  greater  than  2A 
to  break  Cooper  pairs  and  create  excess  quasiparticles. 

To  efficiently  detect  the  excess  quasiparticle  population 
produced  by  photon  depairing  one  must  measure  a  property 
that  is  sensitive  to  all  such  quasiparticles.  These  include  the 
voltage  across  a  series  of  superconductive  weak  links,  the 
microwave  surface  resistance  of  the  superconductor,  and  the 
kinetic  inductance  of  the  superconductor. 

Approaches  pursued  previously  included  the  measure¬ 
ment,  in  granular  HTS  films,  of  the  voltage  signal  induced 
by  laser  radiation,  caused  by  local  suppression  of  the  gap, 
and  thus  the  intergranular  critical  current  in  the  film.  How¬ 
ever,  in  such  a  situation,  only  those  quasiparticles  created 
within  a  distance  of  order  of  the  coherence  length  from  a 
grain  boundary  are  effective  in  reducing  the  intergranular 
critical  current,  and  thus  producing  a  voltage  signal.  The  area 
efficiency  of  such  a  detector  goes  (in  two  dimensions)  as  the 
ratio  of  the  coherence  length  to  the  grain  size  squared 
which  for  a  typical  granular  YBCO  film  might  be 
of  order  10“^.  This  is  an  unacceptably  poor  area  efficiency 
(equivalent  to  a  very  poor  effective  quantum  efficiency)  and 
it  is  incompatible  with  a  high-performance  detector. 

Using  microwave  circuits  we  can  measure  changes  in  the 
microwave  reflectivity*^  and  resonator  Q  produced  with  pho¬ 
toabsorption.  These  detector  approaches  are  sensitive  to  the 
entire  quasiparticle  or  superconducting  condensate  population. 
Such  techniques  are  not  practical  for  detector  array  develop¬ 
ment  because  the  size  of  the  microwave  readout  circuits  is 
incompatible  with  the  small  size  of  detectors,  about  50  X  50 
\Lvr?, 

Another  method  sensitive  to  the  entire  quasiparticle  or 
condensate  population  monitors  the  kinetic  inductance  of  a 
superconducting  film  incorporated  into  a  superconducting 
quantum  interference  device  (SQUID)  loop.  The  period  of 
the  threshold  characteristic  of  gate  current  versus  control 
current  /,.  is  given  by  where  <I>o  is  the  flux  quantum 


and  L  is  the  sum  of  the  geometrical  and  kinetic  inductances 
of  the  SQUID  loop.*  *  Thus,  biasing  at  fixed  values  of  4  and 
4  in  the  voltage  state,  the  voltage  will  be  a  sensitive  function 
of  the  kinetic  inductance,  and  thus  of  the  excess  quasiparticle 
population.  The  disadvantage  of  such  a  scheme  is  that  the 
SQUID  must  operate  in  the  voltage  state,  and  a  calculation 
of  the  sensitivity  of  such  a  detector  based  on  the  measured 
noise  characteristics  of  even  the  best  dc  SQUIDs  yields  ex¬ 
tremely  poor  detectivity. 

A  more  practical  technique,  and  one  that  is  also  sensitive 
to  all  the  quasiparticie’s  produced,  is  to  measure  the  kinetic 
inductance  if  of  the  condensate  in  the  superconductor  under 
irradiation.  The  kinetic  inductance  is  a  measure  of  the  inertia 
of  the  Cooper  pairs  in  a  superconductor  and  is  proportional 
to  the  reciprocal  of  the  Cooper  pair  density  n  (see  Sec.  2). 
With  this  approach,  the  detector  is  maintained  in  the  super¬ 
conducting  current  state  (low  noise)  and  changes  in  its  kinetic 
inductance  are  monitored  as  a  change  in  the  magnetic  field, 
as  detailed  in  Sec.  2. 


2  Quantum  Superconducting  Kinetic  Inductance 
Detectors 

The  description  of  the  quantum  superconducting  kinetic  in¬ 
ductance  detector  (QSKID)  illustrated  in  Fig.  1  is  given  in 
this  section.  We  include  the  basis  for  the  approach  selected, 
the  QSKID  structure,  the  operating  principle,  and  projected 
performance.  Finally,  we  present  an  integrated  QSKID  with 
a  SQUID  readout  circuit. 

2.1  Detector  Approach 

Realizing  a  viable  QSKID  requires  meeting  several  important 
technical  objectives.  First  and  foremost,  the  detector  opera¬ 
tion  and  structure  needs  to  demonstrate  a  quantum  response. 
The  quantum  response  needs  to  occur  with  low  electrical 
noise  realizable  in  a  zero-resistance  superconducting  state. 
Maximum  photoresponse  is  desirable  and  is  obtainable  with 
a  maximum  quasiparticle  lifetime.  The  readout  circuit  of  the 
QSKID  is  to  be  simple,  specifically  a  SQUID  circuit.  The 


Photon 


Fig.  1  QSKID  structure,  represented  here  by  its  equivalent  circuit 
with  magnetic  L  and  kinetic  inductances,  consists  of  a  supercon¬ 
ducting  loop  biases  with  a  dc  current  Iq.  The  bias  current  splits  be¬ 
tween  branch  1  and  2  of  the  superconducting  loop  producing  branch 
currents  and  4,  respectively.  For  identical  branches,  i.e.,  with  the 
same  kinetic  and  magnetic  inductance,  current  is  equal  to  current 
4.  Changing  the  kinetic  inductance  by  photodepairing  in  only  one 
branch  gives  rise  to  a  circulating  photocurrent  A/,  producing  a  de¬ 
tectable  magnetic  field. 
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QSKID  structure  compatible  with  the  aforementioned  ob¬ 
jectives  includes 

1.  no  Josephson  Junctions  or  weak  links  in  the  detector 
structure,  thus  facilitating  manufacturability 

2.  maximum  photoresponse  amplitude 

3.  noise-equivalent  power  (NEP)  superior  to  thermal  de¬ 
tectors  facilitated  by  operation  in  the  zero-resistance 
superconducting  state 

4.  broadband  response,  UV  to  far  infrared  fxm) 

5.  detector  structure  compatible  with  a  100%  fill  factor. 

From  a  manufacturability  point  of  view  it  is  undesirable 
to  use  Josephson  junctions  or  weak-link-based  detectors. 
First,  the  coherence  length  in  high  7^  materials  is  too  short 
for  making  Josephson  junctions,  and  making  reliable  weak 
links  in  high  7^  materials  is  very  difficult.  Second,  the  re¬ 
quired  weak-link  area  for  each  detector  makes  a  large  focal 
plane  impractical.  In  high  T’o  the  coherence  length  ^  is  very 
small  (<2  nm)  and  for  a  photon  to  be  detected  it  needs  to 
be  absorbed  within  a  coherence  length  of  the  junction.  Be¬ 
cause  a  typical  detector  is  about  50  X  50  in  area,  many 
weak  links  would  be  needed  to  provide  a  100%  fill  factor. 
The  weak  links  would  have  to  be  separated  by  about  2  nm 
and  span  the  detector  width.  An  alternative  is  a  vertical  ge¬ 
ometry  where  a  single  weak  link  per  detector  is  used  and  it 
is  oriented  perpendicular  to  the  incident  photons.  Neither 
geometries  are  technologically  practical  if  we  recognize  that 
a  typical  focal  plane  is  expected  to  have  thousands  of  detec¬ 
tors.  The  total  cross-sectional  area  of  the  weak  links  would 
grow  to  greater  than  1  X  1  cm^.  Such  a  detector  array  would 
need  to  satisfy  demanding  critical  current  uniformity  require¬ 
ments.  Junctions  with  a  nonuniform  critical  current  density 
distribution  would  have  significant  fractions  of  areas  with 
zero  current  flow,  because  of  the  small  Josephson  penetration 
depth,  resulting  in  insensitive  regions. 

The  second  issue  deals  with  maximizing  the  photo¬ 
response  amplitude.  It  has  been  mentioned  earlier  that  the 
QSKID’ s  responsivity  is  directly  proportional  to  the  lifetime 
of  the  quasiparticle  excitations.  The  quasiparticle  lifetime  is 
minimized  as  the  temperature  increases  toward  7^-,  thereby 
reducing  the  photoresponsivity.  Quasiparticle  lifetime  re¬ 
duction  occurs  because  as  the  temperature  increases,  so  does 
the  number  of  quasiparticles  and  therefore  the  likelihood  of 
finding  two  with  opposite  momentum  to  recombine  into  a 
Cooper  pair.  Making  use  of  phonon-trapping  phenomena  and 
operating  well  below  Tc,  the  quasiparticle  lifetime  and  thus 
the  photoresponse  amplitude  are  maximized. 

The  third  objective  is  the  use  of  a  detector  structure  not 
operating  in  the  voltage  state.  Satisfying  this  criterion  is  im¬ 
portant  to  achieve  low  operating  power  and  minimum  noise 
operation.  If  we  operate  high  7<^  detectors  in  a  voltage  state 
(with  a  gap  voltage  of  about  20  mV)  and  at  a  bias  current  of 
10  mA,  the  power  consumed  by  each  detector  is  about  0.2 
mW.  With  thousands  of  detectors,  this  becomes  a  significant 
and  an  intolerable  load  on  the  refrigerator.  Furthermore,  given 
the  low  resistance  of  superconductors  operating  in  a  voltage 
state,  the  current  shot  noise  would  be  too  large,  leading  to 
poor  sensitivity.  Using  voltage  readout  circuits,  to  take  ad¬ 
vantage  of  the  small  voltage  noise,  is  unacceptable  because 
of  refrigeration  requirements.  Very  low  noise  voltage  sem¬ 


iconducting  readout  circuits  consume  too  much  power  to  be 
compatible  with  large,  low-power  focal  plane  arrays.  Thus 
if  low  power  and  low  noise  operation  are  sought  the  detector 
needs  to  remain  in  the  zero-resistance  superconducting  state. 

The  fourth  and  fifth  objectives  are  important  from  an  op¬ 
erational  point  of  view.  HTS  materials  have  a  small  energy 
gap  and  therefore  are  expected  to  respond  from  x  rays  to  far 
infrared  (Xc>20  fxm).  The  need  for  a  high  fill  factor  detector 
structure  is  to  simply  maximize  the  detector  quantum  effi¬ 
ciency  and  performance. 

A  QSKID  structure  has  been  formulated'^  to  satisfy  the 
objectives  raised  in  this  section.  The  detector  structure  is 
illustrated  in  Fig.  1 ,  where  a  single  square  superconducting 
loop  forms  the  detector.  No  weak  links  are  used,  and  the 
detector  is  made  simply  of  a  thin  film  patterned  into  a  square 
loop.  The  output  signal  is  a  magnetic  flux  produced  by  the 
circulating  current  in  the  square  loop.  Details  on  the  operation 
of  this  QSKID  are  given  next. 

2.2  QSKID  Structure  and  Operation 

In  this  section,  the  operation  of  the  QSKID  is  described  and 
the  responsivity  is  computed.  First,  the  operating  constraints 
are  computed.  Next,  a  photon  flux  is  used  to  cause  a  pertur¬ 
bation  to  the  quiescent  operating  condition  from  which  the 
responsivity  is  computed. 

The  current  division  between  the  QSKID’ s  two  branches 
is  adjusted  to  minimize  the  total  energy  of  the  detector.  The 
total  energy  (E)  of  the  QSKID  is  made  up  of  the  magnetic 
field  contribution  and  the  kinetic  energy  (£^)  of  the 
supercurrent.  For  the  QSKID  structure  in  Fig.  1,  the  input 
current  Iq  (less  than  the  critical  current)  is  divided  between 
the  two  branches  yielding  currents  /j  and  4.  The  total  mag¬ 
netic  field  energy  Stored  in  the  magnetic  inductances  L, 
and  L2,  of  branches  1  and  2,  respectively,  is 

E^,=\  ■  (5) 

The  kinetic  energy  of  the  supercurrent  4  flowing  in 
branches  1  and  2,  expressed  in  terms  of  the  Cooper  pair 
volume  density  n,  the  Cooper  pair  velocity  v,  the  Cooper  pair 
mass  2m,  the  pair’s  electronic  charge  2e,  the  length  of  each 
branch  €,  and  the  cross  section  of  each  branch  A  is 

£/^=(n|Aj€,)mvf +  (n2^2^2)'”^2  •  (^) 

Recognizing  that  /j=2n,A,evi  and  4  =  2n2^2^'^2’  we  can 
recast  Eq,  (6)  to  obtain  an  expression  for  in  terms  of  /, 
and  4  as 

.  (7) 

The  recasting  of  Eq.  (6)  in  terms  of  supercurrents  /j  and  4 
instead  of  v,  and  V2  results  in  a  form  similar  to  Eq.  (5)  for 
the  magnetic  energy  Em-  The  two  coefficients  iSj  and  ^2 
provide  a  convenient  form  for  expressing  the  kinetic  energy 
of  the  supercurrent  in  terms  of  currents  /,  and  4*  Because 
the  forms  of  Eqs.  (5)  and  (7)  are  similar,  the  coefficients 
and  X2  are  labeled  as  kinetic  inductances  in  correspondence 
to  the  magnetic  inductance  label  given  to  Lj  and  L2  in  Eq. 
(5).  Equations  for  the  kinetic  inductances,  associated  with 
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the  kinetic  energy  of  the  supercurrent  flow  in  a  rectilinear 
geometry  QSKID,  are  defined  as 


\A  jle^n 


The  kinetic  inductance  is  a  function  of  the  detector’s  ge¬ 
ometry  and  is  inversely  proportional  to  the  Cooper  pair  den¬ 
sity  n.  Combining  Eqs.  (5)  and  (7),  the  QSKID’s  total  energy 
E  becomes 


Division  of  current  Iq  —  Ix  +4  between  branches  1  and  2  (in 
Fig.  1)  into  /j  and  is  according  to  the  minimum  energy 
constraint,  or 

bE  dE 

—  =  —  =  (L,+^,)I^-iL2  +  ^2)^2  =  0■  (10) 

Ol I  0/2 

The  QSKID’s  minimum  energy  criterion  and  its  quiescent 
operating  condition  require  the  current  division  between  the 
two  branches  to  satisfy 

+  .  (li) 

A  photoresponse  occurs  in  the  QSKID  when  one  branch, 
branch  1,  for  example,  is  exposed  to  a  photon  irradiance  Ep 
sufficiently  energetic  to  break  Cooper  pairs.  As  Cooper  pairs 
are  broken,  the  kinetic  inductance  if,  will  increase  [see  Eq. 
(8)].  A  change  in  the  kinetic  inductance  of  8if ,  in  branch  1 
changes  by  A/,  the  current  flowing  in  branch  1,  or  equiva¬ 
lently,  produces  a  circulating  current  A/,  in  the  QSKID  loop. 
Taking  the  differential  of  Eq.  (11),  an  expression  for  the 
photoresponse  current  is  obtained.  The  induced  current  A/, 
is  expressed  in  terms  of  the  changes  in  the  kinetic  inductance 
8^,  of  branch  1  as 


AT  =  —  I - - - 

*  *  (Lj +Z.2“b'^2) 


cannot  be  made  equal  to  zero.  Because  the  superconductor's 
wave  function  i|i(A)  is  single  valued,  a  superconducting  cur¬ 
rent  flowing  in  a  closed  loop  F  (see  Fig.  1)  has  a  phase  9 
that  can  only  change  by  ItiN,  where  N  is  an  integer.  Per¬ 
forming  the  closed  loop  line  integration  on  Eq.  ( 1 3),  we  obtain 

jB-dS-hj  =  .  (14) 

.V 

where  B  is  the  magnetic  field  inside  the  detector’s  square 
loop  and  S  is  the  area  enclosed  by  the  curve  F.  Completing 
the  integration  in  Eq.  (14)  and  making  use  of  Eq.  (8),  we 
obtain  the  expression 

Lx/x-L^/2  +  ^x/x-^2^2  =  ^Oo  •  (15) 

Taking  the  differential  /,  in  Eq,  (15)  relative  to  the  changes 
in  the  kinetic  inductance  8if ,  produces  the  same  result  as  given 
by  Eq.  (12),  derived  from  minimum  energy  requirements.  The 
magnetic  field  Bqskid  produced  by  the  circulating  current  A/ 
is  a  measure  of  the  absorbed  photon  irradiance  Ep, 

2.3  QSKID  Responsivity 

The  QSKID’s  responsivity  in  webers  per  watt  is  expressed 
in  terms  of  the  incident  photon  irradiance  Ep.  The  change  in 
the  kinetic  inductance  dif ,  caused  by  this  irradiance  is 


where  n,,  Tj,  T^ff,  and  d  are,  respectively,  the  Cooper  pair 
density,  the  number  of  Cooper  pairs  broken  by  each  photon, 
the  effective  quasiparticle  lifetime,  and  the  detector  thickness. 
The  optical  power  P  incident  on  the  detector  with  an  optical 
area  is  simply  the  product  of  three  terms: 

P=Ap,Ephv  .  (17) 

The  QSKID’s  responsivity  91  in  webers  per  watt  is  computed 
by  combining  Eqs.  (12),  (16),  and  (17)  to  yield 


Note  that  the  photoresponse  current  A/,  can  also  be  computed 
from  fluxoid  conservation  requirements,  and  the  results  are 
the  same  as  given  by  Eq.  (12),  derived  from  minimum  energy 
requirements.  Specifically,  in  a  superconductor  the  conden¬ 
sate  described  by  a  wave  function  that  in  general  can  be 
expressed  by  a  complex  quantity  il;(x)  =  Vn(x)  At  any 
given  point  x  inside  the  superconductor  the  density  of  Cooper 
pairs  is  given  by  n(x),  where  n(x)  =  i|;*(x)i|/(x).  Writing  out 
explicitly  the  current  density  in  a  superconductor  using  the 
canonical  momentum,  we  obtain 


_  2it(  ^  m 


where  d>o  ~  hfle,  A  is  the  vector  potential,  and  is  the  su¬ 
percurrent  density.  The  closed  loop  of  the  detector  is  made 
from  a  superconducting  film  thinner  then  the  magnetic  pen¬ 
etration  depth.  Thus  the  supercurrent  J,  does  not  vanish  any¬ 
where  inside  the  QSKID  loop,  F,  and  the  term  in  Eq.  (13) 


9t(webers/watt)  —  (L,  +  Ly)—r  = 

dP 


^1 


L,+ie,+L2+ig2 


The  responsivity  is  expressed  as  a  product  of  four  factors. 
The  first  factor  is  simply  the  bias  current  /,  flowing  through 
the  branch  illuminated  by  the  photon  irradiance.  This  is  ex¬ 
actly  given  by  Eq.  (11),  and  for  a  balanced  detector  it  is  equal 
to  IqH.  The  second  factor  in  Eq.  (18)  is  that  the  ratio  of 
inductances  is  equal  to  about  unity,  because  the  kinetic  in¬ 
ductances  are  smaller  in  value  than  the  magnetic  inductances. 
The  third  factor  in  Eq.  ( 1 8)  is  a  ratio  of  the  number  of  Cooper 
pairs  broken  by  a  photon  divided  by  the  total  number  of 
Cooper  pairs  in  the  detector.  For  a  0.1-eV  phonon  and  a 
YBCTO  QSKID  50  X  50  p.m  in  area  50  nm  thick  with  a  Cooper 
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pair  density  of  7  X  10“^^  this  term  can  be  estimated 

equal  to  10*'.  Combining  all  these  approximations,  the 
expression  for  the  responsivity  reduces  to 

^li(webers/vvatt)  =  0-5  X  10' ,  (19) 

It  is  evident  that  the  responsivity  is  directly  proportional  to 
the  QSKID  bias  current  /q,  the  effective  quasiparticle  lifetime 
Tgf,,  and  the  kinetic  inductance  ifj.  Substituting  reasonable 
values  for  4==^  I  A  and  nH,  the  responsivity  simplifies 
to  Tgff.  Because  Iq  is  below  the  superconductors  crit¬ 

ical  current,  no  power  is  dissipated  by  the  QSKID.  The  es¬ 
timated  quasiparticle  lifetime  (see  Sec.  2.4)  is  greater  than 
10  fxs  at  10  K  considering  phonon  trapping  effects  in  the 
film.  For  these  values,  the  responsivity  conservatively  com¬ 
putes  to  about  10"''  WbAV  at  10  K.  This  can  readily  be 
detected  by  a  SQUID  circuit. 

2.4  Quasiparticle  Lifetime 

The  photoresponse of  YBCO  and  Nb  were  measured  with 
the  transient  photoimpedance  response  (TPR)  method'"  to 
observe  the  nonequilibrium  photoresponse of  supercon¬ 
ductors,  and  in  particular  to  estimate  the  duration  of  the  non¬ 
equilibrium  excitation’s  lifetime.  Laser  pulses  of  300  fs  at 
665  nm  were  used  to  study  nonequilibrium  excitation  life¬ 
times  from  the  photoinduced  impedance  transient  that  ap¬ 
pears  as  a  voltage  signal  across  the  films.  From  the  amplitude 
and  temporal  dependence  of  the  TPR  signal,  the  thermali- 
zation  process  was  studied  in  superconductors  in  the  normal, 
transition,  and  superconducting  states.  The  temporal  depen¬ 
dence  of  the  thermalization  process  forms  the  basis  for  our 
rough  estimation  of  quasiparticle  lifetime  and  the  viability 
of  the  QSKID. 

The  TPR  photoresponse  clearly  demonstrated  that  the  du¬ 
ration  of  the  photoinduced  nonequilibrium  was  dominated 
by  extrinsic  quasiparticle  lifetime  effects.  Specifically, 
phonon  trapping  in  the  thin-film  samples  significantly  in¬ 
creased  the  quasiparticle  recombination  lifetime  over  intrin¬ 
sic  values.  This  was  particularly  evident  from  TRP  mea¬ 
surement  in  the  normal  state  where  the  signals  were 
bolometric  in  origin  and  phonon  trapping  governed  the  pho¬ 
toresponse  lifetime.  Phonon  trapping  in  thin  films  (<80  nm) 
resulted  in  thermalization  lifetimes  much  greater  than  several 
nanoseconds.  These  lifetimes  were  longer  than  anticipated 
from  a  simple  acoustic  boundary  mismatch  model  by  at  least 
an  order  of  magnitude.  A  nonideal  boundary  between  sub¬ 
strate  and  thin  film  and  the  slower  velocity  expected  for  en¬ 
ergetic  Debye  phonons  in  YBCO  help  to  explain  the  observed 
long  thermalization  lifetimes.  It  is  important  to  stress  that,  in 
the  superconducting  state,  the  quasiparticle  lifetime  will  be 
also  governed  by  phonon  trapping,  and  we  infer  from  the 
data  that  their  lifetime  will  be  much  larger  than  10  ns.  Clearly 
the  net  recombination  of  quasiparticles  cannot  occur  faster 
than  the  escape  time  for  phonons  from  the  film  into  the  sub¬ 
strate  or  the  anharmonic  phonon  decay  lifetime.  The  normal 
state  measurements  put  a  lower  limit  on  the  phonon  escape 
time.  However  in  the  superconducting  state  the  phonon  trap¬ 
ping  time  is  further  lengthened  by  quasiparticle  excitations. 

In  the  superconducting  state,  the  excess  energy  being  ther¬ 
mal  ized  is  shuttled  back  and  forth  between  creating  phonons 
and  exciting  quasiparticles.  This  excess  energy  can  be  re¬ 
moved  from  the  film  only  by  phonon  emission  into  the  sub- 
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strate  or  anharmonic  phonon  decay.  The  energy  stored  in 
quasiparticles  cannot  be  removed  from  the  film  via  quasi¬ 
particle  emission  into  the  substrate  because  of  their  different 
band  gap  structures.  Thus  during  the  time  the  excess  energy 
is  stored  in  quasiparticle  excitations  no  excess  energy  escapes 
into  the  substrate.  At  lower  temperatures,  typically  the  quas¬ 
iparticle  generation  lifetime  T;;  is  much  shorter  than  the  quas¬ 
iparticle  recombination  lifetime  Therefore,  the  excess  en¬ 
ergy  is,  most  of  the  time,  in  the  form  of  excited  quasiparticles. 
Because  the  excess  energy  can  escape  into  the  substrate  only 
through  phonons,  the  storage  of  this  excess  energy  in  the 
form  of  quasiparticles  effectively  increases  the  phonon  trap¬ 
ping  time.  According  to  Rothwarf  and  Taylor,'^  the  effective 
lifetime  for  quasiparticle  recombination  can  be  related  to 
the  phonon  trapping  time  (t.^)  intrinsic  quasiparticle 

generation  and  recombination  lifetimes  as 


"^eff 


(20) 


Because  T;j>Tg,  the  effective  recombination  time  for  quas¬ 
iparticles  should  be  longer  than  the  thermalization  time  con¬ 
stant  for  the  normal  state,  see  Eq.  (20).  From  our  mea¬ 
surement  of  the  normal  and  superconducting  state  responses, 
it  is  evident  that  quasiparticle  lifetimes  much  longer  than  10 
ns  are  likely.  Calculations  indicate  that  the  quasiparticle  life¬ 
times  at  10  K  may  in  fact  be  much  greater  than  10  jis. 

2.5  QSKID  Readout  Circuit 

The  QSKID  readout  circuit  senses  the  magnetic  flux  produced 
by  the  photoinduced  current  circulating  in  the  square  loop 
(see  Fig.  2).  For  achieving  the  highest  sensitivity,  a  SQUID 
is  used  to  sense  the  photoinduced  magnetic  flux.  In  the  de¬ 
tector  currently  under  development,  the  SQUID  is  made  up 
of  two  niobium  tunnel  junctions.  The  QSKID  structure  is 
equally  compatible  with  a  SQUID  readout  circuit  made  from 
YBCO  weak  links.  For  best  sensitivity,  the  SQUID  readout 
circuit  should  be  well  coupled  to  the  flux  produced  by  the 
QSKID,  and  this  is  achieved  by  integrating  the  QSKID  and 
SQUID  on  a  single  substrate.  The  responsivity  of  the  QSKID 
is  about  10“"  WbAV  or  about  0.5  10"  4>o/W,  a  quantity 
readily  detected  by  a  SQUID. 

3  Conclusions 

For  operation  at  cryogenic  temperatures,  excellent  semicon¬ 
ductor  quantum  detectors  are  available.  Intrinsically,  quan¬ 
tum  detectors  are  more  sensitive  than  bolometers.  Thus,  if 
superconducting  detectors  are  to  displace  existing  semicon¬ 
ductor-based  quantum  detectors,  they  need  to  offer  improved 
performance,  possible  only  by  developing  QSKID  and  not 
superconducting  bolometers.  Such  a  superconducting  quan¬ 
tum  detector  has  been  introduced  and  its  operation  described. 
The  projected  responsivity  is  calculated  for  a  SQUID  readout 
circuit.  The  QSKID  is  a  broadband  detector  with  a  response 
extending  beyond  X^20  |jLm  and  consumes  very  little  power. 
These  performance  features  are  not  available  with  semicon¬ 
ducting  detectors.  Ongoing  developments  on  the  theoretical 
and  experimental  front  should  provide,  in  the  near  future, 
answers  to  the  ultimate  performance  of  the  QSKID.  The  in¬ 
formation  presented  represents  a  starting  point  in  the  devel¬ 
opment  of  a  new  superconducting  detector  structure. 
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Fig,  2  QSKID  with  a  SQUID  readout  circuit.  The  QSKID  detector, 
represented  by  the  bottom  square  loop,  is  overlayed  by  a  light  shield 
and  an  insulator.  The  topmost  layer  corresponds  to  the  SQUID  read¬ 
out  circuit,  drawn  enlarged.  The  actual  size  of  the  SQUID  readout 
circuit  Is  much  smaller  than  the  size  of  the  QSKID  detector. 
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The  discovery  of  high-temperature  superconductors  has  led  to  great  efforts  to  find  potential 
applications,  including  the  development  of  photon  detectors.  We  review  the  limitations  of  the 
different  approaches  proposed  for  realizing  superconducting  photon  detectors.  The  structure 
and  operation  of  a  new  quantum  superconducting  kinetic  inductance  detector  (QSKID)  with 
a  SQUID  readout  circuit  is  described.  The  QSKID  is  made  from  a  superconducting  loop 
where  the  photosignals  are  generated  in  response  to  photoinduced  changes  in  the  condensate’s 
kinetic  inductance.  The  QSKID  operates  in  the  zero-resistance  superconducting  state,  thereby 
circumventing  Johnson  noise. 


KEY  WORDS:  HTS;  photodetector;  quantum  detector;  kinetic  inductance. 


1.  INTRODUCTION 

The  discovery  of  high-temperature  superconduc¬ 
tivity  (HTS)  by  Bednorz  and  Muller  [1]  in  1986 
evoked  proposals  for  making  IR  photon  detectors  in 
HTS.  These  proposals  were  offered  because  IR  detec¬ 
tors,  like  HTS,  operate  at  liquid-nitrogen  tempera¬ 
tures.  Enomoto  and  Murakami  [2]  made 
photoresponse  measurements  on  granular 
BaPbo.7Bio.3O3  and  reported  encouraging  results. 
Soon,  more  papers  [3-5]  followed,  reporting  on  the 
photoresponse  of  different  HTS  detector  structures. 
After  much  analysis,  it  was  recognized  [6]  that  the 
reported  photoresponse  signals  were  consistent  with 
a  bolometric  process.  This  was  discouraging  because 
the  sensitivity  [7,8]  of  bolometric  detectors  is  less  than 
the  sensitivity  of  quantum  detectors.  However, 
because  of  technological  limitations  in  semiconduct¬ 
ing  quantum  detectors,  the  search  for  a  superior 
superconducting  quantum  detector  continues.  In  this 
paper,  we  report  on  the  limitations  of  different  pub¬ 
lished  approaches  toward  realizing  HTS  supercon¬ 
ducting  detectors  and  an  approach  to  realizing  a  HTS 
quantum  detector. 


'Westinghouse  Advanced  Technology  Laboratories,  Winterson 
and  Nursery  Roads,  Linthicum,  Maryland  21090. 

^Westinghouse  Science  and  Technology  Center,  1310  Beulah  Road, 
Pittsburgh,  Pennsylvania  15235-5098. 


2.  COlWENTIONAL  SUPERCONDUCTING 
DETECTORS 

Semiconducting  and  superconducting  detectors 
can  be  grouped  into  two  main  categories:  quantum 
and  bolometric  (or  thermal)  detectors.  The  detectors 
realized  in  LTS  and  HTS  materials  can  be  classified 
into  the  aforementioned  categories.  Other  exotic 
mechanisms  of  detection,  such  as  photon-induced 
vortex-antivortex  nucleation,  and  photon-enhanced 
flux  flow,  have  been  proposed  [9]  but  are  not  well 
understood  vis-a-vis  the  four  conventional  detector 
approaches  enumerated  below. 

L  Thermal  detectors,  where  the  photoab- 
sorbed  energy  changes  the  lattice  tempera¬ 
ture  and  some  electronic  property  of  the 
material,  such  as  resistance.  LTS  and  HTS 
transition  edge  bolometers  have  a  sensitivity 
limited  by  the  energy  fluctuations  of  the  lat¬ 
tice  temperature  and  the  photon  flux.  These 
detectors  are  attractive  for  very  long  wave¬ 
length  (>20/im)  applications  where  semi¬ 
conductor-based  detectors  require  cooling  to 
liquid-helium  temperatures. 

II.  Josephson-effect  based  detectors  are  particu¬ 
larly  useful  for  millimeter  and  submillimeter 
wave  detection.  These  devices  use  the 
nonlinear  interaction  between  the  incident 
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photons  and  the  high-frequency  oscillating 
currents  and  voltage,  in  a  biased  Josephson 
junction,  to  perform  millimeter  wave  mixing 
and  direct  detection. 

III.  Quasiparticle  tunneling  devices  such  as  S- 
I-S  mixers  and  direct  detectors,  where  the 
incident  photons  are  absorbed  by  quasipart¬ 
icles,  allowing  tunneling  events  which  would 
otherwise  be  energy-forbidden  to  affect  the 
junction  I-V  characteristics. 

IV.  Pair-breaking  detectors,  where  incident  pho¬ 
tons  (with  energy  larger  than  the  energy  gap) 
break  Cooper  pairs,  leading  to  changes  in 
such  macroscopic  parameters  as  the  super¬ 
conducting  gap  and  the  kinetic  inductance. 
In  these  detectors  photoabsorption  is  by 
quantum  excitations. 

The  basic  requirement  for  a  superconducting 
quantum  detector  are:  (1)  quantum  excitations  to 
create  excess  quasiparticles  by  breaking  Cooper  pairs 
with  photons  that  have  an  energy  greater  than  2A, 
and  (2)  efficient  detection  of  the  excess  quasiparticles 
by  measuring  a  property  sensitive  to  the  condensate 
population.  Properties  sensitive  to  the  quasiparticle 
(or  condensate)  population  include  the  voltage  across 
a  series  of  weak  links,  the  superconductor’s  micro- 
wave  surface  resistance,  or  the  superconductor’s 
kinetic  inductance. 

In  granular  HTS  films,  the  voltage  signals  are 
induced  by  laser  radiation  that  locally  suppresses  the 
gap  and  the  intergranular  critical  current.  Only  those 
quasiparticles  created  within  a  distance  of  the  order 
of  the  coherence  length  ^  from  the  grain  boundaries 
produce  a  voltage  signal.  The  area  efficiency  of  a  gran¬ 
ular  HTS  detector  is  expressed  (in  2D)  as  the  ratio  of 
the  coherence  length  to  the  average  grain  size,  Lg, 
squared.  For  a  typical  granular  YBCO  film,  this  area 
efficiency  ratio  {^jLg'f  is  of  the  order  10”^.  This  is 
an  extremely  poor  area  efficiency  and  equivalent  to  a 
very  poor  quantum  efficiency,  incompatible  with  high- 
performance  detectors. 

The  poor  area  efficiency  in  granular  HTS  detec¬ 
tors  can  be  overcome  with  microwave  readout  circuits 
which  are  sensitive  to  the  entire  quasiparticle  or  con¬ 
densate  population,  such  as  changes  in  the  microwave 
reflectivity  [10]  and  resonator  Q  produced  with 
photoabsorption.  The  microwave  readout  circuit’s 
size  is  much  larger  than  a  typical  pixel  size,  about 
50x50^m^.  Thus,  such  a  readout  approach  is  not 
practical  nor  compatible  with  the  development  of 
dense  two-dimensional  arrays. 


Another  method  sensitive  to  photoinduced 
changes  in  the  condensate  population  is  the  voltage 
vs.  flux  signal  of  the  superconducting  film  incorpor¬ 
ated  into  a  SQUID  loop.  Biasing  this  SQUID  loop 
in  the  voltage  state,  photodepairing  changes  the  con¬ 
densate’s  kinetic  inductance,  resulting  in  a  voltage  sig¬ 
nal,  derived  from  the  SQUID’s  periodicity  with  flux. 
The  period  of  the  SQUID’s  threshold  characteristic 
with  applied  flux  is  given  by  where  ^>0  is  the 

flux  quantum  and  L  is  the  total  inductance,  equal  to 
the  sum  of  the  geometrical  and  kinetic  inductances  in 
the  SQUID  loop  [11].  Such  a  detector  must  operate  in 
the  voltage  state,  and  calculations  reveal  a  sensitivity 
severely  limited  by  the  SQUID’s  noise  characteristics, 
resulting  in  poor  detectivity. 

To  realize  sensitive  (low-noise)  detectors,  it  is 
critical  that  the  LTS  or  HTS  detectors  operate  in  the 
current  and  not  voltage  state.  This  detector  needs  to 
be  sensitive  to  photoinduced  changes  in  the  con¬ 
densate,  e.g.,  sensitive  to  the  condensate’s  kinetic 
inductance  if.  With  such  an  approach,  the  detector 
remains  in  the  superconducting  current  state  (low 
noise),  and  changes  in  its  kinetic  inductance  are  moni¬ 
tored  as  a  change  in  the  magnetic  field,  as  detailed  in 
the  next  section. 

3.  QUANTUM  SUPERCONDUCTING  KINETIC 
INDUCTANCE  DETECTORS 

A  viable  QSKID  needs  to  operate  in  the  super¬ 
conducting  current  state  to  achieve  low-noise  per¬ 
formance.  For  maximum  photoresponse,  long 
quasiparticle  photoexcitation  lifetimes  are  necessary. 
The  QSKID’s  read-out  circuit  must  be  simple, 
specifically  a  SQUID  circuit.  The  QSKID  structure 
compatible  with  the  aforementioned  objectives 
includes :  (a)  no  weak  links,  facilitating  low  noise  and 
manufacturability,  (b)  maximum  photoresponse  and 
NEP  by  operation  in  the  zero-resistance  supercon¬ 
ducting  state,  (c)  multispectral  response  from  UV  to 
the  far  infrared  (Ac>20  /im),  and  (d)  structure  com¬ 
patible  with  a  100%  fill  factor.  Meeting  these  objec¬ 
tives  is  necessary  to  produce  a  viable  QSKID  with 
a  performance  superior  to  competing  semiconductor 
quantum  detectors.  Such  a  detector  structure  has  been 
formulated  [12]  and  is  schematically  illustrated  in  Fig. 

1  by  a  single  superconducting  loop.  No  weak  links 
are  used,  and  the  output  signal  is  the  magnetic  flux 
produced  by  the  circulating  current  in  the  loop. 
Details  on  the  QSKID’s  operation  are  given  next. 

The  division  of  current  h  <  h  into  currents  /j  and 
h  satisfies  an  energy  minimum  [12,13]  and  the  single¬ 
valuedness  for  the  condensate  wave  function. 
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Fig.  1.  The  QSKID  (a)  and  its  equivalent  circuit  (b)  consists  of 
superconducting  loop  biased  with  a  dc  current  /o,  with  magnetic 
(L)  and  kinetic  (J5f)  inductances.  Bias  current  /o  splits  between 
branches  1  (A )  and  2  (/2)  of  the  superconducting  loop.  For  ident¬ 
ical  branches  h-h.  Photodepairing  only  in  branch  1  changes  the 
kinetic  inductance  {5^\)  and  gives  rise  to  a  circulating  photo- 
current  A/,  producing  a  detectable  magnetic  field.  A  small  SQUID 
placed  (not  shown)  over  the  opening  in  the  QSKID  detects  the 
photoinduced  magnetic  flux. 


expressed  by  a  complex  quantity  \f/(x)  =  y/n(x)  We 
calculate  the  current  division  by  requiring  a  single¬ 
valued  if/ix),  and  the  result  is  consistent  with  energy 
minimization  requirements.  At  any  point  x  inside  the 
superconductor,  the  density  of  Cooper  pairs  is  given 
by  rt(x),  where  n{x)  =  y/*(x)y/(x).  Writing  out  explic¬ 
itly  the  condensate’s  phase  gradient  V^,  we  obtain 


^  Ztt  ^  m  - 

<I>o  \  2e^n 


area  {a)  is 


m 

le^n 


(4) 


In  a  QSKID  with  a  symmetrical  geometry  shielded 
from  external  fields  and  not  subjected  to  any  illumina¬ 
tion,  the  net  enclosed  flux  is  zero  [the  right  side  of 
Eq.  (3)  is  zero]  and  no  current  circulates,  /i=/2  = 
Iq/1,  Illuminating  only  branch  number  1  in  Fig.  1 
changes  The  change  in  that  will  produce 
a  circulating  current  A/,  obtained  by  taking  the 
differential  of  Eq.  (3),  is 


-  (5) 

The  circulating  current  A/  produces  a  magnetic 
field  ^qskid  and  flux  0  =  A/(Li+L2),  which  is  a 
measure  of  the  absorbed  photon  flux  <I>ph.  The 
QSKID’s  responsivity  9?  (Webers/watt)  (dependent 
on  the  flux  Oph,  changes  in  if  i ,  and  number  of  broken 
pairs  )  is 


5if,= 


d0ph/  «i 


(6) 


where  wi,  t?,  rcfr,  and  d  are,  respectively,  the  Cooper 
pair  density,  the  number  of  pairs  broken  by  a  photon, 
the  effective  quasiparticle  lifetime,  and  the  detector 
thickness.  The  incident  power  P  on  the  detector,  with 
an  optical  area  Ax,  is 


where,  <^o  =  h/2e,  A  is  the  vector  potential,  and  Js  is 
the  supercurrent  density.  In  superconducting  films 
thinner  than  the  London  penetration  depth,  the  super¬ 
current  Js  does  not  vanish  anywhere  inside  the 
QSKID  loop,  and  \i/{x)  is  single-valued.  The  line  inte¬ 
gral  in  Eq.  (1),  after  some  simplifications,  yields 


P^Ax<^p^.hv  (7) 

Combining  Eqs.  (5)  and  (6),  multiplying  by  (Lj  +  L2), 
and  dividing  by  7,  the  computed  QSKID’s  responsiv¬ 
ity  9?  (Webers/watt)  becomes 

9?  (Webers/watt) 


B  dS+  t  -^Js  dl=N<bo  (2) 

J  ^  2r/i 

where  B  is  the  magnetic  field  inside  the  detector’s  loop 
and  S  is  the  area  enclosed  by  the  integration  curve. 
Completing  the  integration  in  Eq.  (2),  we  obtain 

Lxlx  - LJi  +  if  1/.  -  JP2I2  =  N^o  (3) 

where  the  four  coefficients  if|(Li)  and  if 2(^2)  are 
the  kinetic  (magnetic)  inductances  of  branch  1  and  2, 
respectively.  The  kinetic  inductance  of  each  branch 
expressed  in  terms  of  the  Cooper  pair’s  mass  (2m), 
charge  {2e),  the  branch  length  (/),  with  cross  sectional 


_ Y  T]  y.g^iTefl-'j 

L\+^s+L2  +  ^2)\Axn\d)\  hv  } 

The  first  factor  is  the  bias  current  I\  flowing 
through  the  illuminated  branch,  and  for  a  balanced 
detector  it  is  equal  to  h/l.  The  second  factor,  the 
ratio  of  inductances,  is  equal  to  about  unity,  because 
and  £f2~L2.  The  third  factor,  for  a  0.1-eV 
photon  and  a  YBCO  QSKID  50  fjm  x  50  tim  in  area 
and  160nm  thick  with  a  Cooper  pair  density  of 
«10^‘’cm“\  is  estimated  to  be  2  x  10“".  Combining 
all  these,  the  expression  for  the  responsivity  reduces 
to 


91  (Webers/watt)  s  lO^/oTerri^’i . 


(9) 
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The  responsivity  is  proportional  to  the  QSKID 
bias  current  /q,  the  effective  quasiparticle  lifetime  icfr, 
and  the  kinetic  inductance  if  j .  Substituting  values 
for  a  and  %2  nH,  the  responsivity  simpli- 

fies  to  9t^0.02reflr.  Since  Iq  is  below  the  superconduc^ 
tor’s  critical  current,  no  power  is  dissipated  by  the 
QSKID.  The  estimated  quasiparticle  lifetime  [14,15] 
is  about  10  msec  at  10  K,  considering  phonon  trap¬ 
ping  effects  in  the  film.  For  these  values  the  responsiv¬ 
ity  conservatively  computes  to  about  lO""^  Webers/ 
Watt  at  10  K.  This  can  readily  be  detected  by  a 
SQUID  circuit  with  a  sensitivity  of  %10“^’  Webber/ 

The  QSKID  read-out  circuit  senses  the  magnetic 
flux  produced  by  the  photoinduced  current  circulating 
in  the  loop  (see  Fig.  1).  A  SQUID  is  used  to  sense 
the  photoinduced  magnetic  flux.  In  the  QSKID  under 
development,  the  SQUID  is  made  up  of  two  niobium 
tunnel  junctions  or  YBCO  weak  links.  For  best  sensi¬ 
tivity,  the  SQUID  read-out  circuit  should  be  well 
coupled  to  the  flux  produced  by  the  QSKID,  integrat¬ 
ing  the  QSKID  and  SQUID  on  a  single  substrate. 
The  responsivity  of  the  QSKID  is  about  10“"^  Webers/ 
Watt  or  about  0.5  x  lO'*  Oo/Watt,  a  quantity  readily 
detected  by  a  SQUID,  yielding  an  NEP^^r^lO"^’ 
Watt/Vffe. 

4.  CONCLUSIONS 

Excellent  semiconductor  quantum  detectors  are 
available  for  operation  at  cryogenic  temperatures. 
Intrinsically,  quantum  detectors  are  more  sensitive 
than  bolometers.  Thus,  if  superconducting  detectors 
are  to  replace  existing  semiconducting  quantum  detec¬ 
tors,  they  need  to  offer  improved  performance,  pos¬ 
sible  only  by  developing  QSKID  and  not 
superconducting  bolometers.  A  superconducting 
quantum  detector,  with  a  SQUID  read-out  circuit,  is 


introduced  and  its  operation  is  described.  The  QSKID 
projected  responsivity  is  calculated.  The  QSKID  is  a 
multispectral  detector  with  a  response  extending 
beyond  A^20/im  and  consumes  very  low  power. 

ACKNOWLEDGMENTS 

This  research  was  supported  by  Westinghouse 
IR&D  under  task  Z6I601EHAA,  the  Naval  Research 
Laboratory  under  tasks  No.  N0014-92-C-2033,  and 
the  Ballistic  Missile  Defence  Organization/Innovative 
Science  &  Technology  Office  under  task  No. 
DASG60-92-0140, 

REFERENCES 

1.  J.  G.  Bednorz  and  K.  A.  Muller,  Z.  Phys.  B  64,  189  (1986); 
Rev.  Mod.  Phys.  60,  585  (1988). 

2.  Y.  Enomoto  and  T.  Murakami,  J.  Appl.  Phys.  59,  3807  (1966). 

3.  H.  S.  Kwok,  J.  P.  Zheng,  and  Q.  Y.  Ying,  Appi.  Phys.  Lett. 
54,  2473  (1989). 

4.  M.  Leung,  P.  R.  Broussard,  J.  H.  Claassen,  M.  Osofsky,  S.  A. 
Wolf,  and  U,  Strum,  Appl.  pAys.  Lett.  51,  2046  (1987). 

5.  W.  S.  Brocklesby,  D.  Monroe,  A.  F.  J.  Levi,  M.  Hong,  S.  H, 
Liou,  J.  Kwo,  C.  E.  Rice,  P.  M.  Mankiewich,  and  R.  E. 
Howard,  Appl.  Phys.  Lett.  54,  1175  (1989). 

6.  M.  G.  Forrester,  M.  Gottlieb,  J.  R.  Gavaler,  and  A.  I.  Bragin- 
ski,  IEEE  Trans.  Magn.  MAG-25,  1327  (1989). 

7.  R.  A.  Smith,  F.  E.  Jones,  and  R.  P.  Chasmar,  The  Detection 
and  Measurement  of  Infra-red  Radiation  (Oxford  University 
Press,  New  York,  1958),  p.  206. 

8-  P-  W.  Kruse,  L.  D,  McGIauchlin,  and  R.  B.  McQuistan,  Ele¬ 
ments  of  Infrared  Technology:  Generation,  Transmission,  and 
Detection  (J,  Wiley,  New  York,  1962),  p.  354. 

9.  A.  M.  Kadin,  M.  Leung,  A,  D.  Smith,  and  J.  M.  Murduck, 
IEEE  Trans.  Magn.  27,  1540  (1991). 

10.  R.  Kaplan,  W.  E.  Carlos,  and  E.  J.  Cukauskas,  J.  Appl.  Phys. 
67,  4212  (1990). 

11.  E.  K.  Track,  M.  Radparvar,  and  S.  M.  Paris,  IEEE  Trans. 
Magn.  MAG-25,  1096  (1980). 

12.  N.  Bluzer,  United  States  Patents  5,179,072  and  5,185,527. 

13.  N.  Bluzer  and  M.  G.  Forrester,  Proceedings  of  the  SPIE  confer¬ 
ence,  San  Diego,  California,  July  11-16,  1993. 

14.  N.  Bluzer,  Phys.  Rev.  B  44,  10222  (1991);  N.  Bluzer,  Phys. 

46,  1033  (1992). 

15.  N.  Bluzer,  /.  Appl.  Phys.  71,  1336  (1992);  IEEE  Trans.  Appl. 
Supereond.  3,  2869  (1993). 


Presented  at  the  1994  ASC  meeting  in  Boston,  MA  and  published  in: 
IEEE  Trans,  on  Applied  Superconductivity,  5(2),  2583  (1995). 


Quantum  Detectors  in  Superconducting  YBCO* 

Nathan  Bluzer 

Westinghouse  Advanced  Technology  Labs.,  Winterson  &  Nursery  Rds.,  Linthicum,  MD  21090  USA 

Martin  G.  Forrester 

Westinghouse  Science  and  Technology  Center,  1310  Beulah  Rd.,  Pittsburgh,  PA  15235  USA 


Abstract — A  superconducting  quantum  detector  structure  in 
YBCO  is  described  with  a  directly  coupled  SQUID  read  out 
circuit.  The  detector  geometry  is  optimized  for  maximum 
photoresponse  with  the  use  of  a  serpentine  pattern.  The 
serpentine  pattern  reduces  quasiparticle  diffusion  effects  and 
maximizes  the  photoinduced  changes  in  the  kinetic  inductance. 
The  operation  of  this  sensor  is  analyzed  in  terms  of  geometry 
and  quasiparticle  lifetime  to  obtain  expressions  for  the  signal 
and  noise  of  this  detector.  The  background  limited  infrared 
performance  of  this  detector  is  calculated  to  be  about  5x10’*^ 
NEP.  The  photoresponse  is  expected  to  be  limited  by  the 
superconducting  energy  gap,  about  30-40|im. 

1.  Introduction 

Since  the  discovery  of  High  Temperature  Supercon¬ 
ductivity  (HTS)  by  Bednorz  and  Muller  [1]  many 
proposals  were  offered  for  fabricating  Infrared  (IR) 
detectors  in  HTS.  Enomoto  and  Murakami  [2]  measured 
the  photoresponse  in  BaPbo^Bio  3O3  and  reported 
encouraging  results.  Many  more  photoresponse 
measurements  were  reported  on  HTS,  however  it  was 
shown  by  Forrester  et  al.  [3]  that  the  results  reported  were 
consistent  with  a  bolometric  and  not  a  quantum  (or 
photon)  detector  response.  This  was  disappointing  for  it 
has  been  long  recognized  that  theoretically  above  IK 
quantum  detectors  offer  superior  sensitivity  over 
bolometric  detectors  [4].  Hence  the  expectation  of 
replacing  semiconductor  based  quantum  detectors  by 
bolometric  superconductor  based  detectors  was 
unreasonable.  If  superconducting  detectors  are  to 
compete  with  semiconducting  quantum  detectors, 
superconducting  detector  structures  are  needed  which 
provide  photon  noise  limited  performance. 

At  Westinghouse  we  have  been  pursuing,  for  several 
years,  the  development  of  superconducting  quantum 
detectors.  The  superconducting  quantum  detectors  are 
intended  for  imaging  applications.  In  such  applications 
the  detector  integration  times  are  from  several  psec  to 
several  msec  long.  And  thus  we  alleviate  the  need  for 
very  fast  readout  circuits  and  very  fast  photoresponse 
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measurements  to  short  laser  pulses.  The  detectors  are 
intended  to  operate  in  the  superconducting  current  state 
(zero  voltage)  under  static  nonequilibrium  conditions  and 
in  the  Meissner  state  (below  Hd),  to  circumvent  noise 
from  fluxoid  formations  and  movement.  The  detectors’ 
performance  is  dependent  on  the  quasiparticle  lifetimes 
and  these  lifetimes  are  maximized  as  the  detectors’ 
operating  temperature  is  lowered  significantly  below  the 
superconducting  transition  temperature  Tq.  Additionally, 
the  approach  selected  [5]  does  not  use  weak-links  in  the 
detector  structure  and  therby  maximizes  the  quantum 
efficiency,  and  reduce  the  fabrication  difficulties 
associated  with  making  weak-links  in  YBCO.  Weak-link 
based  detectors  are  only  sensitive  to  photoexcitations  that 
occur  within  several  coherence  lengths  ^  from  the  weak- 
link.  Thus  for  good  quantum  efficiency  ,  many  weak- 
links  need  within  each  detector  pixel  because  ^  «1.5  nm 
is  many  times  smaller  than  the  typical  detector  pixel,  25- 
100pm.  The  use  of  many  weak-links  for  each  detector 
imposes  severe  fabrication  yield  problems  as  well  as 
sensitivity  issues.  Hence  we  have  selected  a  simpler 
detector  and  readout  approach  for  realizing 
superconducting  quantum  sensors. 

II.  Superconducting  Sensor 

The  superconducting  quantum  detector  we  are 
developing  are  based  on  the  condensate’s  kinetic 
inductance  and  the  changes  in  this  kinetic  inductance 
produced  by  photodepairing  of  Cooper  pairs.  Generically, 
the  structure  of  Quantum  Superconducting  Kinetic 
Inductance  Detector  (QSKID)  is  shown  in  Fig.  1,  where 
a  thin  (less  than  a  London  penetration  depth) 
superconducting  closed  loop,  operating  in  the  current  and 
Meissner  states,  divides  a  DC  current  Iq  =  Il+Ir  into  a  left 
branch  current  II  and  a  right  branch  current  Ir.  The 
current  division  satisfies  the  Hamiltonian’s  minimum 
energy  conditions  [5]  [6]  and  is  expressed  by  Equation  1 
as: 

LlIl  =  (1) 

Where  the  L^  and  Lr  are  the  total  inductance  terms  for 
the  left  and  right  branches,  respectively.  The  total 
inductance  of  a  branch  includes:  the  geometrical  (or 
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magnetic  field)  self  inductance,  the  interbranch  mutual 
inductance,  and  the  kinetic  inductance.  For  example,  the 
total  left  branches  inductance  Ll=Llg  -L^m+L^k  includes 
the  geometrical  self  inductance  the  mutual 

inductance  between  the  left  and  right  branches  (with 
negative  coupling)  Ll^,  and  the  kinetic  inductance  Llk- 


Figure  1.  The  generic  diagram  of  a  QSKJD  consisting  of  a 
superconducting  loop  dividing  a  DC  current  Iq  into  and  Ir  .  Only 
one  branch,  the  left,  e.g.,  is  exposed  to  a  photon  flux  while  the  other  is 
shielded.  Changes  in  the  kinetic  inductance  of  the  exposed  branch 
causes  a  circulating  current  51. 

A  photon  flux  'F  on  the  left  branch  breaks  up  Cooper 
pairs  and  produces  a  decrease  in  the  Cooper  pair  density 
Pl  .The  kinetic  inductance  Llk^EIlA^cpI/ESl  PlI  is 
inversely  dependent  on  Pl,  the  length  of  the  left  branch  1^, 
the  mass  of  a  Cooper  pair  mcp,  the  cross  sectional  area  of 
the  left  branch  Sl,  and  the  Cooper  pair  charge  qcp 
squared.  With  photoabsorption  the  Cooper  pair  density  in 
the  exposed  branch  decreases  and  the  kinetic  inductance 
increases  by  ALlk  The  signal  81  developed  with 
photoabsorption  in  the  left  branch  is  calculated  by  taking 
the  differential  of  Eq.  1  to  obtain: 


5/  = 


^LK  ^0 


(2) 


The  photoinduced  signal  81  depends  on  the  QSKID 
geometry  and  quasiparticle  lifetime  and  these  are 
addressed  below. 


A.  Detector  Structure 

With  a  serpentine  geometry  shown  in  Fig.  2  the 
performance  of  the  QSKID  is  improved  over  the  simple 
superconducting  loop  shown  in  Fig.  1.  The  improvements 
occur  because  for  a  given  photon  flux  'F  the  change  in 
ALlk  is  larger  and  because  improved  isolation  is  achieved 
between  the  illuminated  and  masked  branches 


A  serpentine  geometry  QSKID  with  N  vertical 
segments  in  each  branch,  has  the  lengths  of  the  left  and 
right  branchs  IgL  and  Isr  increased  about  N  times  over  a 
nonserpentine  QSKID.  For  the  QSKID  in  Fig.  2,  IsJ  1l«7 
and  IsR  /1r«7.  A  similar  scaling  occurs  for  the.  left  and 
right  cross  sectional  areas  Ssl  and  Ssr  resulting  in 
Ssl/Sl«Ssr/Sr«7.  Also  the  bias  current  Iq  decrease  by  a 
factor  of  7,  to  provide  a  desirable  reduction  for  operating 
integrated  circuits.  The  serpentine  geometry  results  in  a 
significant  increase  in  the  kinetic  inductance  relative  to 
the  geometrical  inductance.  Combining  all  these  factors, 
we  calculate  that  for  the  serpentine  QSKID  under  a  given 
photon  flux  'F,  the  signal  81  increases  approximately  7 
fold. 


S(^  detector  bias  input 


Figure  2.  The  YBCO  superconducting  quantum  detector  is  shaped  into 
a  serpentine  pattern  with  the  directly  coupled  SQUID  readout  circuit 
located  at  the  bottom.  The  photon  flux  only  illuminates  the  left  half  of 
the  serpentine  patterned  detector. 

An  additional  benefit  with  the  serpentine  structure  is 
interbranch  isolation  to  maintain  maximum  signal.  With 
photoillumination  Llk  f^e  quasiparticle  population 
increase.  Interbranch  quasiparticle  and  Cooper  pair 
diffusion  between  the  illuminated  and  the  masked 
branches  will  reduce  the  signal.  Interbranch  diffusion  is 
facilitated  with  quasiparticle  having  long  lifetimes  and 
this  occur  by  phonon  trapping  and  low  temperature 
operation.  For  YBCO,  the  quasiparticle  diffusion  length 
is  estimated  to  be:  2pm  for  Insec,  10pm  for  lOnsec,  and 
60pm  for  Ipsec  lifetimes  respectively.  These  diffusion 
lengths  are  significant  relative  to  the  QSKID  pixel  size.  A 
serpentine  structure  significantly  increase  the  interbranch 


SQUID  flux  null  output  current  node 
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path  lengths  and  thereby  resulting  in  improved 
interbranch  isolation,  i.e.  less  signal  loss. 

B.  Detector  *s  Photoresponse 


Photoabsorption  in  the  left  branch  of  the  serpentine 
QSKID  in  Fig.  2  produces  a  change  in  the  inductance 
The  photoresponse  is  calculated  in  terms  of  the  change  in 
Ljl  using  the  electrical  equivalent  circuit  in  Fig.  3.  The 
equivalent  circuit  includes  two  loops  each  supplied  with  a 
DC  current  generator.  The  top  loop,  supplied  by  current 
lo,  represents  the  serpentine  detector  structure.  Below  it,  a 
second  loop,  supplied  by  a  DC  current  IgQ,  represents  the 
detector’s  SQUID  readout  circuit.  Each  inductor  in  the 
equivalent  circuit  includes:  the  self,  the  mutual,  and  the 
kinetic  inductance  terms  associated  with  the  structure  in 
Fig.  2.  The  equivalent  circuit  is  symmetrical  and  the  left 
(right)  side  inductors  and  weak  links  are  indicated  by  a 
subscript  which  includes  the  letter  L  (R). 


'R 


Figure  3.  Equivalent  circuit  diagram  for  a  directly  coupled  SQUID 
readout  circuit  to  a  YBCO  superconducting  quantum  detector.  There 
are  two  DC  bias  currents,  one  for  the  detector  Iq,  and  one  for  the 
SQUID  readout  device  Isq. 


The  photoresponse  is  calculated  by  relating  the 
photoinduced  currents  in  the  detector  51  with  the  SQUID 
current  5i.  Before  illumination,  and  under  initial 
conditions,  currents  51  and  5i  are  set  equal  to  zero. 
Minimum  energy  conditions  [5]  [6],  represented  by  Eq.  1, 
yield  the  equilibrium  conditions  equation  as: 


[Al  h^2L  ”  [A/?  ■*"  A/?  ]  A  A/? 


The  total  SQUID  loop  inductance  (L2L+L2R+L3L+L3R) 
is  designed  to  be  small  to  achieve  good  critical  current 
modulation  with  flux.  For  QSKID  with  a  large  serpentine 
geometry  results  in  Lil»  L2L  and  Lir»L2r  Therefore, 


currents  i^  and  iR  are  shunted  to  ground  primarily  through 
L2L  and  L2R  ,  respectively,  and  we  can  neglect  the  flow  of 
currents  11  and  Ir  through  Lil  and  Lir,  Integrating  around 
the  SQUID  loop  and  requiring  a  zero  or  a  multiple  of  In, 
the  phase  change,  we  obtain  an  equation  for  the  SQUID 
under  equilibrium  conditions,  i.e. 

—  0^  +  Az,4  Az,[A  4]  ” 

”  o.  .  , 

2^  ® R  A/? [A  A] 

Where  Oq  is  a  fluxoid,  0l  and  0r  represent  the  phase 
shift  across  the  left  and  right  weak  links. 

With  the  application  of  a  flux  'F,  photodepairing 
increases  Ljl  by  5L  and  establishes  new  static  conditions 
with  currents  51  and  51.  circulating  in  the  detector  and 
SQUID  loops,  respectively.  Making  the  following 
substitutions  in  Eq.  3:  Il=Il-5I,  Ir=Ir+5I,  Lil=L,l+5L, 
iL=iL+5I,  iR=iR+5i,  and  keeping  only  first  order  terms  ,  the 
nonequilibrium  equation  for  the  detector  is: 


d  LIl  -  5/  [L^j^  +  +  P\L  +  ^2L  ] 


(5) 


Incorporating  the  same  substitutions  into  Eq.  4  and 
keeping  only  the  first  order  terms,  the  nonequilibrium 
equation  for  the  SQUID  becomes: 


^50  +  5/  [^3;^  +  +  Lji  +  Z,2^  ] 

n 

=  6/  [L2R  +  L2i\ 


(6) 


Eq.  5  is  readily  solved  by  recognizing  that  since 
Lil»  L2L  and  Lir»L2r  the  last  term  can  be  neglected, 
and  we  obtain  a  solution  similar  to  Eq.  2.  Similarly,  Eq.  6 
is  solved  by  recognizing  that  when  operating  in  the 
voltage  state,  the  leading  term  is  the  smallest  since  on 
average  50«O  and  the  phase  angle  across  the  weak  link 
remains  very  close  to  nil.  Combining  Eq.  5  and  6  we 
obtain  a  relationship  between  the  SQUID  and  detector 
signal  currents,  specifically, 

6/ =5/ - -  (7) 

[^R  +  A  +  Ax  +  Ai?] 

The  signal  coupling  efficiency  between  the  QSKID’s 
51  and  the  SQUID’s  5i  is  optimized  by  making  the 
relative  sizes  of  [L3R+L3L]  «  [L2L+L2R].  The  detector 
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geometry,  illustrated  in  Fig.  2,  is  designed  for  achieving  a 
maximum  coupling  efficiency 

The  squid’s  voltage  signal,  produced  by  the 
QSKID’s  signal  51,  is  obtained  from  the  RSJ  model  [7] 
[8]  in  the  limit  as  the  weak  link’s  shunting  capacitance  C 
goes  to  zero.  The  SQUID ’s  average  output  voltage  signal 
5V,  with  weak  links  that  each  have  an  equivalent 
shunting  resistance  of  R,  is: 

=  - g;  (8) 

Va^-1 

where  a=IsQ/2Ic  and  Ic  is  equal  to  the  weak  link’s  critical 
current.  The  sensitivity  of  the  QSKID  is  examined  next. 

C  Detector 's  Sensitivity 

Since  the  mechanism  of  HTS  is  not  well  understood, 
it  is  difficult  to  provide  accurate  predictions  of  static 
nonequilibrium  phenomena  in  these  materials.  However, 
using  the  BCS  model  and  the  Rothwarf-Taylor  equations, 
sensitivity  predictions  have  been  made.  An  extensive 
analysis  has  been  made  [9]  [10]  and  here  we  shall  only 
give  highlights  of  these  results 

The  noise  in  the  QSKID,  due  to  fluctuations  in  the 
Cooper  pair  population,  is  calculated  to  be  much  less  than 
the  Johnson  noise  in  voltage  state  SQUID  The  QSKID 
sensitivity  is  limited  by:  the  current  noise  in  the  SQUID 
readout  circuit  This  noise  can  be  minimized  by  using 
SQUIDs  that  have  a  smaller  Ic  and  a  larger  loop 
inductance,  but  at  the  cost  of  noise  rounding  in  the 
SQUID’S  I-V  characteristics. 

The  photoresponse  was  modeled  for  Background 
(photon  noise)  Limited  Infrared  Performance  (BLIP) 
where  the  photon  noise  is  larger  than  the  detector  and 
readout  circuit  noises,  added  in  quadrature.  Hence,  the 
ultimate  performance  level  is  achieved  by  minimizing  the 
detector  and  the  readout  circuit  noises.  For  BLIP,  the 
photogenerated  quasiparticle  population  needs  to  be 
larger  than  the  thermal  equilibrium  quasiparticle 
population.  Also,  the  detectors  responsivity,  dependent  on 
the  quasiparticle  lifetime,  needs  to  be  maximized  to 
overcome  readout  circuit  noise.  The  51  is  maximized  by 
operating  the  detector  at  lower  temperatures  (30K  e.g.) 
and  utilizing  phonon  trapping  to  realize  large  effective 
quasiparticle  lifetimes.  In  YBCO  Pair  Breaking  (PB) 
phonon  are  trapped  because  anisotropy  naturally  confines 
the  movement  of  PB  phonons  in  the  ab  plane.  Such 
phonon  trapping  is  limited  by  anharmonic  decay  lifetime 
of  PB  phonons  into  non-pair  breaking  phonons. 


We  have  calculated  the  BLIP  sensitivity  of  a 
lOOxlOOjLim  QSKID  with  a  directly  coupled  SQUID 
readout  circuit.  The  calculated  sensitivity,  for  a  lOOpm 
square  detector  operating  in  a  30Hz  bandwidth,  in  terms 
of  Noise  Equivalent  Power  is  about  10’^^- 10’^^  Watts. 
This  calculated  level  of  performance  is  limited  by  the 
detector’s  SQUID  readout  circuit.  The  QSKID  projected 
spectral  response  is  limited  by  the  Cooper  pair  binding 
energy,  predicted  for  YBCO  to  be  about  30-|im 

III.  Conclusions 

A  QSKID  structure  with  a  directly  coupled  SQUID 
readout  circuit  were  described  and  the  performance 
analyzed.  The  operation  of  this  sensor  structure  is  based 
on  photoinduced  changes  in  the  QSKID  kinetic 
inductance.  A  serpentine  pattern  is  used  for  maximizing 
the  photoresponse.  The  detectors  sensitivity  is  limited  by 
the  current  noise  in  the  directly  coupled  SQUID  readout 
circuit.  The  predicted  BLIP  NEP  is  10’^^- 10"^^  Watts  with 
a  cut-off  wavelength  between  30-40pm.  The  predicted 
performance  makes  the  QSKID  attractive  and 
competitive  particularly  for  detection  at  very  long  IR 
wavelength  of  cold  objects. 
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